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Abstract 

A High Spectral Resolution Lidar (HSRL) that uses an iodine absorption filter and 
a tunable, narrow-bandwidth NdrYAG laser is demonstrated. The iodine absorption 
filter provides better performance than the Fabry-Perot etalon that it replaces. T his 
study presents an instrument design that can be used as a basis for a design of simple 
and robust lidar for the measurements of the optical properties of the atmosphere. 
The HSRL provides calibrated measurements of the optical properties of the atmo- 
spheric aerosols. These observations include measurements of aerosol backscatter 
cross section, optical depth, backscatter phase function, depolarization, and multiple 
scattering. The errors in the HSRL data are discussed and the effects of different 
errors on the measured optical parameters are shown. 
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1 Introduction 

The climate models used for modeling the transport of the short and longwave radia- 
tion in the atmosphere require a knowledge of the aerosol and cloud optical properties. 
Because the optical properties of cirrus clouds are not well known, the measurements 
of cloud optical depth, phase function, and particle size provide important informa- 
tion. Clouds affect the radiative balance of the earth and its atmosphere by reflecting 
incoming solar radiation and trapping outgoing longwave radiation. Cirrus clouds 
have been found to have an important role on this process' ,2 . Cirrus clouds consist 
of large, nonspherical ice crystals and they are generally found between altitudes of 
4 and 15 km. Compared to water clouds, cirrus clouds are generally optically thin. 
Studies based on climate models suggest that presence of cirrus may produce either 
a heating or cooling effect depending on the cloud optical properties and altitude 1 ’ 2 . 

The University of Wisconsin High Spectral Resolution Lidar (HSRL) measures 
optical properties of the atmosphere by separating the Doppler-broadened molecu- 
lar backscatter return from the unbroadened aerosol return. The molecular signal is 
then used as a calibration target which is available at each point in the lidar pro- 
file. This calibration allows unambiguous measurements of aerosol scattering cross 
section, optical depth, and backscatter phase function. Also measurements of depo- 
larization and multiple scattering can be performed. In this study, clear air aerosols, 
stratospheric aerosols, and cloud particles are all referred to as aerosols. Similar 
measurements of cloud optical parameters can be made with a Raman-lidar 3,4 , but 
because the Raman backscatter cross section is about 1000 times smaller than the 
Rayleigh backscatter cross section, the HSRL has a significant signal strength ad- 
vantage over the Ram an- technique. Another advantage of the HSRL is that it can 
provide daytime measurements while sky noise background limits the measurements 
of the weak Raman signal to night time. 

The basic idea of an HSRL was originally presented by Fiocco and DeWolf 5 . 
They proposed the measurements of atmospheric aerosols by interferometrically sep- 
arating the backscatter signal. They demonstrated the principle of an HSRL with 
laboratory measurements of scattering from natural aerosols and artificially produced 
dense fog. A later experiment with the HSRL technique was performed by Schwiesow 
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and Lading 6 , who used Michelson interferometers in an attempt to measure atmo- 
spheric temperature. Their evaluation showed that a Michelson interferometer based 
measurements would theoretically produce accurate measurements of atmospheric 
temperature, but development of a functional system was limited by the quality of 
available optical components. 

An investigation done by Shipley et a l. 7 for a shuttle borne lidar experiment to 
measure global distribution of aerosols and their effects on the atmospheric heat 
budget started the University of Wisconsin HSRL research. The demonstration of 
the first University of Wisconsin HSRL was published by Shipley et a i. 8 . The paper 
by Sroga et a l. 9 presented the first results measured with the same system. The 
transmitter was based on a dye laser that operated at 476.8 nm wavelength. This 
laser provided only a 2-4 mW output and it operated at 100 Hz repetition rate. The 
backscatter signal was collected by a 0.35 telescope with a 320 /rrad receiver field of 
view. The filtering of the solar background was performed by using an interference 
filter and three low resolution etalons with a total bandpass of 2 pm. The separation 
between aerosol and molecular backscatter signals was based on a high resolution 
etalon with ~ 0.5 pm bandpass. The HSRL was mounted in an aircraft. The system 
pointed down to ground and it therefore covered only about 2 km measurement range. 
The range dependence of the receiver spectrometer bandpass due to the angular 
sensitivity of the etalon transmission complicated the measurements. 

The aircraft system was later modified for ground based operations and the 
changes were reported by Grund 10 . The problems due to the range dependence of 
the spectrometer bandpass were reduced by using a fiber optics scrambler 12 . The 
transmitter laser was changed to a CuCl laser operating at 510.6 nm wavelength. 
The laser output power was 50 mW at 8 kHz repetition rate. The receiver telescope 
and the high resolution etalon remained the same, but the number of low resolution 
etalons was decreased to two providing a 2.5 pm bandpass. This system was ca- 
pable of measuring cirrus cloud optical parameters up to ~ 12 km altitude. Later, 
the development of lasers enabled the use of an injection seeded, frequency doubled 
Nd:YAG- laser (532 nm) as the HSRL transmitter. The new laser provided a 0.7 W 
output at 4 KHz repetition rate. This change was reported by Grund and Eloranta 11 . 
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After all the development in the HSRL, problems with the insufficient calibration 
accuracy and environmental sensitivity of the system operation remained. The HSRL 
produced accurate measurements of clear air and thin cirrus, but the measurement 
accuracy was not sufficient for the cases where the aerosol signal was large compared 
to the molecular signal. Eventhough the measurement accuracy was not high enough 
to provide direct measurements of the optical depth profile inside a cloud, the total 
cloud optical depth was obtained. The cloud optical depth was determined from 
the decrease in the molecular signal across the cloud. The calculations of optical 
depth profile inside a cloud were made by assuming a constant backscatter phase 
function with altitude. The measured molecular profile was fitted with a profile 
derived from the known molecular backscatter cross section that had been corrected 
for the atmospheric extinction. The extinction profile was obtained by applying the 
Bernoulli solution to the aerosol backscatter profile 12 . The measured optical depth 
was used as a constraint. 

The accuracy of the HSRL calibration was improved and the sensitivity for envi- 
ronmental changes was reduced in the next University of Wisconsin HSRL (Eloranta 
and Piironen 13 ). In order to make the system operation more stable, the system was 
moved to a temperature controlled and vibrationally isolated environment. Partici- 
pation in field experiments was made possible by moving the lidar to a semitrailer. 
The diameter of the receiver telescope was increased to a 0.5 m and the field of view 
was decreased to 160 /rrad. A new photon counting data system was designed and ex- 
tensive computer control of the system operations was implemented. Depolarization 
and multiple scattering measurement capabilities were added. With the improved 
system, the measurements of cloud optical properties were performed at altitudes up 
to 30 km. 

During the 18 first years of the University of Wisconsin HSRL, the separation be- 
tween aerosol and molecular scattering was based on a Fabry-Perot etalon with a 0.5 
pm bandpass. With this etalon, the signal was divided into two channels. The signal 
reflected from the etalon surface was measured with one channel. This channel con- 
tained the the wings of the molecular spectrum and the part of the aerosol spectrum 
that did not pass the high resolution etalon. The other channel was used to measure 



4 


the part of the aerosol spectrum and the central part of the molecular spectrum that 
transmitted through the etalon. The basic idea of the spectral separation of a high 
resolution etalon based HSRL is presented in Fig 1. 


Etalon Transmittance Transmitted Spectrum 



Figure 1. The spectral response of a high resolution etalon based HSRL. The received 
signal is a combination of the Doppler-broadened molecular backscatter spectrum and the 
unbroadened aerosol spectrum. The aerosol signal and the center of the molecular signal 
that tr ansmi ts through the etalon is detected with one channel. The other channel detects 
the signed reflected from the etalon. 

For the aerosol backscatter signal, only a ^ 2:1 separation ratio between aerosol 
and molecular channel is achieved by a etalon. The Fabry-Perot etalon based HSRL 
produces accurate measurements of clear air, thin cirrus, and stratospheric aerosols. 
However, when the system is used to probe dense water clouds, the aerosol signal 
becomes on the order of 10^ larger than the molecular return. Therefore, the inversion 
coefficients used to separate the aerosol and molecular signals must be known with 
better than 0.01% accuracy or otherwise some of the aerosol return will appear in 
the separated molecular return. Since the inversion coefficients for the etalon based 
system are known with only ~ 0.1% accuracy, the measurements of dense water 
clouds are subject to error. The performance of a Fabry-Perot etalon is limited by its 
finesse and the angular distribution of incoming light. The etalon must be operated 
in pressure and temperature controlled environment, since better than a 0.1 mbar 
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pressure tuning accuracy and 0.1 °C temperature stability are required to keep the 
filter performance stable. 

Shimizu et al. 14 proposed the use of an narrow- band atomic absorption filter in an 
HSRL and She et al. 15 reported high spectral resolution lidar measurements of tem- 
perature and aerosol extinction coefficient made by using a barium atomic absorption 
filter. These studies, and later studies from the same group 16,17 , have shown, that 
an absorption filter provides a high rejection against aerosol scattering and therefore 
it makes the separation between molecular and aerosol scattering easier. Another 
advantage of an absorption filter is the stability of the absorption characteristics 14 . 
Furthermore, the transmission characteristics of an absorption filter are not depen- 
dent on the mechanical alignment of the filter 14 or the angular dependence of the 
incoming light. Also, a wide dynamic range in rejection against aerosol scattering is 
achieved by simply changing the vapor pressure 14 or the length of the cell. 

This study presents an HSRL employing an iodine absorption filter. The spectrum 
of the £ 3 n+ -4 A 1 !^ electronic transition in molecular iodine has more than 22 000 
absorption lines in the visible wavelengths 18 , and 8 of them are easily reached by 
thermally tuning a frequency doubled Nd:YAG laser output 19 . Compared to the 
barium, the advantage of iodine is that instead of requiring a dye laser, a narrow 
bandwidth, frequency doubled Nd:YAG laser can be used. Also, a strong absorption 
is obtained in a short cell at room temperature. Even though iodine has extensive 
hyperfine structure, the absorption line width is similar to the barium line width, 
which is broadened by operating at a temperature of ~500 °C. 

A large number of iodine absorption lines have been used as reference for Doppler- 
limited spectroscopy studies and also numerous spectroscopic studies of the line struc- 
ture and hyperfine structure have been performed. Liao and Gupta 20 reported a use of 
an iodine absorption vapor cell as a narrow band filter for fluorescence spectroscopy. 
Lately Miles et al. 21 reported the measurements of flow field properties based on an 
iodine absorption filter and Filtered Rayleigh Scattering technique. The first iodine 
absorption filter based HSRL is presented here. 

The basic idea of an iodine absorption filter based HSRL is presented in Figure 
2. The received backscatter signal is divided to two channels. One channel detecting 
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a sample from the total backscatter spectrum and the other channel the spectrum 
filtered by the iodine absorption filter. This signal contains information about the 
wings of the molecular spectrum and a small aerosol cross-talk signal. The first mea- 
surements made with the iodine absorption based HSRL were presented by Piironen 
and Eloranta 22 . 


Transmitted Spectrum 



Figure 2. The spectral response of an iodine absorption filter based HSRL. The detected 
backscatter spectrum is divided to two channels. One containing the information from the 
total backscatter spectrum and the other the wings of the molecular spectrum and a small 
aerosol cross-talk. 

The work presented in this thesis is organized as follows. A short introduction to 
the HSRL theory is given in Chapter 2. After that, the instrumentation of the new 
iodine absorption based HSRL is given in Chapter 3. This chapter also presents the 
basic principles of the polarization and multiple scattering measurements. A more 
detailed description of the iodine absorption cell is given in Chapter 4. The system 
calibration and laser wavelength locking are discussed in Chapter 5. In Chapter 6, 
examples from the data obtained with the new HSRL are given. This chapter gives 
a more detailed description of the depolarization measurements and also shows the 
effects of multiple scattering on the measured depolarization ratio. In addition to the 
depolarization measurements, the measurements of scattering ratio, aerosol backscat- 
ter cross section, and optical depth are discussed. The measurements show, that the 
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use of an iodine absorption filter enables accurate measurements of cloud optical pa- 
rameters. Because the cross talk between channels can be accurately corrected and 
because the 160 /irad field of view of the HSRL effectively suppresses multiple scat- 
tering, the optical depth inside a cloud can be measured. This makes future studies of 
scattering phase function possible. As a final example from the HSRL measurements 
the Chapter 6 shows an atmospheric temperature profile obtained by the HSRL. The 
accuracy of the HSRL measurements is discussed in Chapter 7 and the error analysis 
is presented. 
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2 Lidar theory of the HSRL 

The lidar equation of a monostatic lidar system can be presented in the following 
form 


= fr + A/(fl) + B, 


(i) 


where 

N(R) 

N. 

c 

A 

0a(R),Pm(R) 

P g (?T ,/?) Prn (tt ,/£) 

4tt 5 47 r 

t(R) 


M(R) 

B 


= number of photons incident on the receiver field of 
view from range R per data bin length, 

= number of transmitted photons, 

= speed of light, 

= area of the receiver telescope, 

= aerosol and molecular scattering cross sections per 
unit volume, respectively, 

= normalized backscatter phase function due to aerosol and 
molecular scattering, respectively, 

= one way optical depth between lidar and backscatter 
volume at range R, 

= fo* A(r')dr', where 

0t(R) = total extinction cross section per unit volume, 

= multiple scattering contribution incident on the receiver 
per data bin length, 

= number of background photons incident on the receiver 
per data bin length. 


Equation (1) demonstrates that the lidar return depends on both the local value of 
the backscatter cross section and on the optical depth between the lidar and backscat- 
tering volume at range R. A conventional single channel lidar provides a single mea- 
surement of N(R) at each range, and therefore it does not provide enough information 
to solve the lidar equation for extinction or backscatter cross section. The knowledge 
about the relationship between the backscatter cross section and extinction has to 
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be known or assumed and an estimate of a boundary value has to be given (see for 
example Klett 23 ). 

In order to make measurements of extinction and backscatter cross section, the 
HSRL measures two signals which can be processed to present lidar returns from 
aerosol and molecule backscattering. The separation between the aerosol and molec- 
ular scattering is possible, because the signal backscattered from air molecules is 
Doppler-broadened, while the signal from more massive, slowly moving aerosols re- 
mains spectrally unbroadened. 

In the case of the HSRL, two separate lidar equations can be written 

1) For aerosols 

= F(R)^-HM) p ‘^ - e~ 2 f" ( 2 ) 

2) For molecules 

N m (R) = F(R)^t3 m (R)±e-*l>- (3) 

The molecular backscatter phase function is here replaced with its analytical value 

The term F(R) is an overlap term, that is a function of receiver and transmitter 
geometry 8 . Because the molecular scattering cross section is a function of the air 
density and can be calculated from the Rayleigh scattering theory using an indepen- 
dently measured atmospheric density profile, the equation (3) is well defined and it 
can be solved for the extinction. The molecular scattering cross section provides a 
calibration target which is available at each point of the lidar return. 

Following optical parameters can be obtained from the separated lidar signals: 

1) Aerosol to molecular backscatter ratio (scattering ratio) SR(R) 


SR(R) = 


Naim 

N m (R)’ 


(4) 


Notice that, the backscatter ratio measurement is not dependent from the overlap 
term, and therefore measurements of backscatter ratio for the overlap region can be 
performed. The same effect can be seen later on the aerosol backscatter cross section 
and depolarization ratio determination. 
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2) Aerosol backscatter cross section 

The properties of the molecular scattering are well known and therefore, the at- 
mospheric density can be used to compute the molecular scattering cross section per 
unit volume 0 m . Prom Rayleigh scattering theory 24 , the molecular scattering cross 
section per unit volume is 


where 

M m 

d n 


(5) 

= number of gas molecules per unit volume 

= differential Rayleigh scattering cross section at scattering angle i r 
per average gas molecule 


For mixture of atmospheric gases below 100 km altitude 24 , 


doRjjr) 

dQ 


5.45[ 


A(pm)._, 


0.55 


pio 


-28 


2_— 1 


cm sr 


(6) 


The number of gas molecules per unit volume can be calculated from the atmo- 
spheric pressure P(R) and temperature T(R). 


where 

Kf ( E>\ _ 

Mm{R) ~ T(R ) ft. 

P(R) 

= atmospheric pressure at range R 

T(R) 

= atmospheric temperature at range R 

Na 

= Avogadro’s number 

R a 

= gas constant of dry air 


Thus, by using the previous equations, the molecular scattering cross section per 
unit volume can be written in the following form. 


C air P(R) 

T(R) 


where 


( 8 ) 
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C'air = 3.786 x 10 ~ s KhPa~ i m~ l at 532 nm (from Eq. 5) 

P(R) = atmospheric pressure 

T(R ) = atmospheric temperature 

The aerosol backscatter cross section can be calculated by using the backscatter 
ratio and the calculated value of /3 m 


= S R(R)P m (R)± 


( 9 ) 


3) The optical depth of a layer between ranges R x and R 2 


-(R 2 )-t(Rx)= [ R2 p ( (r')dr' = kn 


/3 m (j7 2 )^lV m (/7 1 )F(fi 1 ) 


0m(Ri)R%N m (R 2 )F(R 2 )\ 


( 10 ) 


4) Total extinction cross section 

The average value of total extinction cross section is a range derivative of the 
optical depth 


dr(R) t(R 2 )-t(R 1 ) 
dR R 2 - Rx 


( 11 ) 


5) Backscatter phase function 

The backscatter phase function, — a ^ R ^ , can be calculated from the following equa- 
tion 


where 


Raj*, R) 

4n 


a(R) = /3 m (R) 


3 SR(R) 
87 r 0 e (R) ’ 


a(R ) = single scattering albedo. 


( 12 ) 


For water and ice clouds and most aerosols, the single scattering albedo at the 
wavelength of 532 nm can be assumed to be unity and therefore the backscatter phase 
function can be measured. 
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6) Linear depolarization ratio 

The discrimination between water and ice clouds can be made from the depolar- 
ization properties of the backscatter signal. The linear depolarization ratio can be 
given as follows. 


$a,m{R) — N(a,m)±(R)/N( a ,m)\\(R), (13) 

where N ± and N { \ are the number of incident on the receiver perpendicular and 
parallel to the receiver polarization axis, respectively. In the case of the HSRL, 
separate depolarizations of aerosol and molecular backscattering can be measured. 

7) Cloud particle size 

By measuring the signal strength variations as a function of field of view, the size of 
the cloud particles can be measured. The multiple scattered lidar return is a function 
of receiver field of view, particle size, range from lidar and the optical depth of the 
cloud. In the HSRL, the multiple scattering information along with the simultaneous 
single scattering measurement are used to calculate the cloud particle size.A detailed 
description of the multiple scattering approximation used for the HSRL measurements 
is presented by Eloranta 25 ' 26 and Eloranta and Shipley 27 . 



13 


3 Instrumentation 

3.1 Introduction 

3.1.1 An overview of the high resolution etalon based system 

In the earlier versions of the University of Wisconsin HSRL, the separation between 
aerosol and molecular backscatter signals was made with a high resolution etalon with 
~ 0.5 pm bandpass. The Figure 3 shows the receiver of the etalon based system 13 . The 
backscattered light was collected by the receiver telescope and directed through a fiber 
optics scrambler. The fiber optics scrambler was used to reduce the range dependence 
of the receiver spectrometer bandpass due to the angular sensitivity of the etalon 
transmission 12 . After collimation, the signal was prefiltered with an interference filter 
and a pair of low resolution etalons. After passing the dual aperture, the light was 
directed to the high resolution etalon. The high resolution etalon was slightly tilted 
with respect to the optical axis. This allowed the light that did not pass through 
the etalon to be reflected back to the dual aperture and to the molecular channel 
photodetector (PMT1). The light that passed through the high resolution etalon was 
detected with the aerosol channel photodetector (PMT2). The signal detected with 
the aerosol channel was a combination of aerosol backscatter spectrum and the center 
of molecular backscatter spectrum. The signal detected with the molecular channel 
was a combination of the wings of the molecular spectrum and the part of the aerosol 
backscatter spectrum that did not pass the high resolution etalon. 
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Figure 3. The receiver of the high resolution etalon based system. The backscatter signal is 
collected with a telescope. The fiber optics scrambler reduces the range dependence of the 
receiver spectrometer bandpass due to the angular sensitivity of the etalon transmission. 
The signal is background filtered with a interference filter and a pair of low resolution 
etalons. The separation between aerosol and molecular backscatter signals is done in the 
high resolution etalon. The part of the spectrum that transmits the etalon is combination 
of aerosol backscatter spectrum and the center of molecular backscatter spectrum (PMT2). 
The signal detected with PMT1 is a combination of the wings of the molecular spectrum and 
the part of the aerosol backscatter spectrum that did not pass the high resolution etalon. 

3.1.2 A short introduction to the improvements 

The instrument has been modified to measure polarization and a separate channel has 
been added to measure both polarization and signal amplitude variations as function 
of receiver field of view. These modifications allow discrimination between ice and 
water clouds and measurements of multiple scattering. Along with these modifica- 
tions, the high resolution etalon, used in the earlier versions of the HSRL, has been 
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replaced with an iodine absorption filter. Also, the system has been redesigned for 
use in a semitrailer. Other changes include new data acquisition system, tempera- 
ture controlled environment for the receiver optics, and improved optical design. The 
system parameters of the HSRL are summarized in Table 1. 

The changes to the HSRL have improved the measurement capabilities of the 
HSRL. With the new HSRL, the measurements of the optical parameters of the 
atmosphere below 35 km can be made and therefore, the measurements can cover 
clouds from water clouds to high altitude cirrus clouds and also measurements of 
stratospheric aerosol layers can be performed. Clouds with optical depths up to ^ 3 
can be measured. This means that most of the cirrus cloud cases can be fully observed 
and the cloud bases of the thick water clouds can be measured up to 300-500 m 
inside the cloud. This has been achieved by using a high pulse repetition rate, small 
pulse energy per laser pulse, and very fast photon counting electronics. Furthermore, 
the averaging time required for a good signal to noise ratio has been decreased by 
improving the optical transmission of the system. The use of the iodine filter, the 
controlling of the pressure of the etalons and the temperature of the optics, and 
the locking of the laser wavelength to an iodine absorption peak have increased the 
stability and the reliability of the system. Therefore, the need for frequent calibrations 
is eliminated and a stable operation can be maintained over a long period of time. The 
enhanced calibration technique has improved the accuracy of the HSRL calibrations. 
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Table 1. HSRL specifications (April 26, 1994) 

Transmitter 


Wavelength 

532 nm 

Pulse duration 

~ 130 ns 

Pulse repetition rate 

4 kHz 

Average power 

~ 0.2 W (seeded) 

Frequency stability: 


Without ^-locking 

~ 100 MHz/hour (= 0.09 pm/h) 

With / 2 -locking 

A± 0.052 pm 


Receiver 


Telescope 
Type 
Diameter 
Focal length 

Wide field of view channel 
Filter bandwidth 
Field of view 
Polarization rejection 

Aerosol+molecular channel 
Filter bandwidth 
Field of view 
Polarization rejection 

Molecular channel 

Filter bandwidth 
Field of view 
Polarization rejection 

Photon counting data system 

Number of counters 
Number of data buffers 
Minimum bin width 
Number of bins 
Max. rate of counters 
Embedded computer 
System control 


Dall-Kirkham 
0.5 m 
5.08 m 

1.0 nm 

computer adjustable 0.21 -4 4.0 mrad 
~lxl0~ 3 


8 pm 

0.160 mrad 
~ lxlO -3 


1.8 pm 
0.160 mrad 
~lxl0 -3 


3 (wfov, aerosol+molec., molec.) 
6 (2 polarizations per channel) 
100 ns 

8192/channel 

~ 10 9 Hz, tested to 350 MHz 
Intel i960CA 
Sun Sparc II 
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3.2 Transmitter 

3.2.1 Transmitter laser 

As a transmitter, the University of Wisconsin HSRL uses a continuously pumped, 
Q-switched, injection seeded, frequency doubled Nd:YAG laser operating at a 4 kHz 
pulse repetition rate 12 (Figure 4). The host cavity is based on a Quantronix model 116 
laser. The quarter waveplates on either side of the Nd:YAG-rod are used to generate 
circularly polarized light in the rod in order to eliminate the spatial hole burning. 
The narrow bandwidth single frequency operation is achieved by injecting the cw- 
output of a diode pumped monolithic Nd:YAG ring laser (a modified Lightwave model 
S-100 seedlaser) into the host cavity through a partially transmitting rear mirror. 
Measurements have shown, that less than 100 MHz/h frequency drift is achieved. 
The 4 kHz pulse repetition rate is achieved by using an acousto-optic Q-switch. 

The output of the laser is externally doubled by a KTP crystal. The frequency 
doubled output of the laser is tunable over a 124 GHz frequency range by controlling 
the temperature of the seedlaser diode. The original analog temperature controlling 
circuit of the seedlaser has been modified and connected to the microprocessor to 
allow remote control of the seedlaser Nd:YAG crystal temperature. 


Cavity ’ump Cavity Focusing 

Mirror Lamp Polarizer Q-Switch Mirror Lens 



Data System Sync. 


Figure 4. The HSRL transmitter laser. 
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After frequency doubling, the 532 nm and 1064 nm beams are separated by a 
harmonic beam separator. Because up to 15% of the residual 1064 nm beam is 
reflected from the beam separator, a second harmonic beam separator in the laser 
output (Figure 5) is used to provide more efficient filtering. The remaining 1064 nm 
beam is cleaned by a spatial filter. The spatial filter also cleans a halo of broader 
bandpass radiation. 

Successful injection seeding requires that overlap between seedlaser resonance fre- 
quency and the frequency of a host cavity longitudinal mode is achieved. When the 
frequency of a host cavity longitudinal mode is locked to the seedlaser resonance 
frequency, the pulse developing out from the seedlaser signal will saturate the homo- 
geneously broadened gain medium preventing development of any other axial modes 
from the spontaneous emission. Therefore, because one longitudinal mode is ampli- 
fied more than any other mode, a spectrally narrow bandwidth pulse is generated. 
Since the seedlaser emissions used for injection seeding are generally ~ 10 orders of 
magnitude stronger than the spontaneous emission, the Q-switched pulse builds up 
sooner out of the seed emission than the spontaneous emission. The frequency lock- 
ing between seedlaser and host cavity is realized by controlling the host cavity length 
with a piezoelectric translator. 

Because the seeded pulse builds up more rapidly when the host cavity is tuned to 
the seed laser frequency, the time between the Q-switch opening and the subsequent 
laser pulse (the Q-switch buildup time) can be used to servo control the tuning of 
the host cavity 28 . The spectral purity of the outgoing laser pulse can also be moni- 
tored using the Q-switch buildup time. A microcontroller based feedback loop seeks to 
maintain the cavity length by dithering the rear mirror so that the minimum Q-switch 
buildup time is obtained. For each laser pulse, a small offset voltage is applied to the 
piezoelectric translator and the effects to the Q-switch buildup time are simultane- 
ously monitored. Based on the Q-switch buildup time the mirror position is driven 
towards the minimum time. Statistics on the Q-switch buildup time is collected at 
4 kHz rate, but the mirror position is dithered at ~ 140 Hz. The amount of dither 
is calculated based on the seeding percentage. The observed Q-switch build up time 
is ~ 4.5 fjs and the difference between seeded and unseeded conditions is typically 
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^ 500 ns. Because the spectral purity of the transmitted laser beam is important for 
HSRL measurements and because some shots are unseeded, the information from the 
Q-switch build up time is used to trigger the data system only for seeded shots. 

The frequency stability of a single- frequency laser is determined by the seedlaser 
and the stability of the frequency locking is affected by the environment. The injection 
seeded operation is easily interrupted by mechanical vibrations transmitted through 
the support structure or through the cooling water or hoses. The effects of change in 
the optical length of the cavity due to thermal expansions of the support structure 
and temperature changes in the laser rod are compensated by the active controlling 
of the main cavity length. The laser is installed on a Super-Invar breadboard sup- 
ported by a thick honey comb table that is mounted into a vibration insulated frame. 
The invar breadboard and the honey comb table are isolated from each other by a 
rubber pad. The laser was delivered with a Super-Invar table, which is thermally 
stable, but it is found to be acoustically sensitive. The active controlling of the main 
cavity length already minimizes the effects of changing environmental temperature 
so that an acoustically more insensitive table might provide a better performance. 
The mechanical vibrations coming from the laser cooling pump and transmitted into 
the cavity by hoses and water are isolated by a pressure dumper in the cooling water 
line. Because the original pressure dumper did not offer enough isolation, an extra 
damping system was installed. Also extra long elastic plastic water hoses are used to 
further attenuate the vibrations coming from the laser cooling pump. Furthermore, 
the temperature of the surrounding environment is stabilized by controlling the air 
temperature around the laser. The wavelength locking of the laser output to a iodine 
absorption line is demonstrated later in Chapter 5.2. 
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Figure 5. The HSRL transmitter setup. The HSRL transmitter employs an injection 
seeded, frequency doubled Nd:YAG laser. A Pockels cell in the output is used to rotate 
the polarization of the outgoing laser beam by 90 degrees for alternative laser pulses. An 
iodine absorption cell is used for the frequency locking of the laser wavelength. A sample 
of each laser pulse is directed to a pair of optical fibers, delayed, and injected back to the 
receiver for system calibrations. The length of the fibers is set so that the time-separated 
pulses can be recorded into the data profile. 


3.2.2 Optics 

A small amount (~2 %) of each transmitted laser pulse is directed into a pair of 100 m 
long optical fibers, delayed, and injected back to the receiver for system calibrations 
(see Figure 5). Since some of the laser shots are unseeded, the delay is necessary 
so that the measured Q-switch buildup time of each pulse can be used as a quality 
control to trigger the data system only for seeded shots. The length of the fibers is 
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set so that the time-separated pulses can be recorded in the data profile. The other 
ends of the calibration fibers are connected to a diffuse cylinder, which is located 
at the receiver such that the calibration light signals can be used to monitor system 
performance during data taking. Another 2 % of the laser light is taken into an energy 
monitor. The energy normalization of the received signal is realized by storing the 
energy of each laser pulse into the data record. A 4 cm long iodine cell is used for 
frequency locking of the laser. This cell provides ~ 50 % absorption when the laser is 
tuned to the absorption peak. The frequency locking to the iodine peak is described 
later in Chapter 5.2. 

Because the HSRL measures depolarization of the clouds and uses this information 
to separate between water and ice clouds, the polarization stage of the outgoing laser 
beam has to be well controlled. The output of the laser is linearly polarized and 
the orientation of the polarization vector is set perpendicular to the plane of the 
hypotenuse of the polarization cube. The use of the polarization cube guarantees 
that the residual cross-polarized component is cleaned out from the outgoing laser 
beam. In order to be able to use the same receiver optics and the same detector for 
both polarization components, the polarization of the outgoing laser pulse is rotated 
by 90 degrees for alternative laser pulses by a Pockels cell. The accuracy of the Pockels 
cell rotation is measured by observing the light coming from the Pockels cell through 
a polarizing cube with a photodetector. These calibrations show that the residual 
cross polarization can be reduced to 0.1 % of the parallel component. The calibration 
accuracy is limited by the 0.1 % rejection accuracy of the polarizing cube. The timing 
of the voltage switching between laser pulses is synchronized with the Q-switch signal 
and the controller electronics assures that the proper Pockels cell voltage is applied 
in time before the next laser shot. The final alignment between the transmitter and 
receiver polarization axes is performed by adjusting a half-waveplate and by using 
the atmosphere as a calibration target. 

The use of small receiver field of views is made possible by decreasing the diver- 
gence of the outgoing laser beam with a beam expanding telescope (magnification 
4x). 
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3.3 Receiver 

3.3.1 Optics 

The backscattered light is collected using a 50 cm diameter telescope (see Figure 
6). To avoid specular reflection from horizontally oriented ice crystals, the system 
pointing direction is tilted by 4 degrees from the zenith. The signal polarized parallel 
to the receiver polarization axis and the signal perpendicular to that are separated by 
using a polarization cube (see Figure 6). This cube separates the signal to the WFOV- 
channel (PMT 3) and to the spectrometer-channels (combined aerosol -(-molecular 
channel (PMT1) and molecular channel (PMT 2)). An extra polarization cube for 
both channels is used to clean up the cross polarization component of the received 
signal down to a level that makes accurate depolarization measurements possible. 

The multiple scattering properties can be studied by using the WFOV-channel, 
while simultaneously measuring the single scattering return with the spectrometer 
channels. A stepper motor driven aperture for the WFOV-channel enables field of 
view variations between 0.224 and 4 mrad. An interference filter with 1 nm bandpass 
is used for the background filtering. Because the field of view can be wide and the 
filter bandpass is broad, the WFOV measurements are currently limited to night time. 
The field of view of the spectrometer channels is 160 /xrad, which further decreases 
the background sky noise and multiple scattering effects to the measured signal. 

In order to achieve daytime measurements with the spectrometer channels, a low 
resolution etalon pair and an 1 nm interference filter are used to suppress the back- 
ground solar radiance. The bandpass of the etalon pair is ~9 pm, which is about 3 
times the width of the expected Doppler broadened molecular spectrum. The length 
of the spacers is chosen so that only one common transmission order for the etalon 
pair occurs inside the interference filter bandpass. The bandpass of the low resolu- 
tion etalon pair is wider than in earlier systems 8,12 in order to decrease the system 
sensitivity for the drift between transmitter wavelength and the etalon transmission 
maximum. The effects of environmental changes to the etalon performance are elim- 
inated by mounting each etalon into a separate pressure and temperature controlled 
chamber. A third etalon chamber is used as a pressure reference. The absolute pres- 
sure in reference chamber and the differential pressures in each etalon chamber are 
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maintained by a computer controlled servo loop. The etalons are operated near the 
pressure where the changes in tuning gas density with temperature are cancelled by 
the thermal expansion of the etalon spacers 29 . Nitrogen is used as a tuning gas. 


To / From Atmosphere 



Calibration 
Light Source 


Field Stop A/4 


Figure 6. A polarizing prism at the output of the telescope separates the orthogonally 
polarized signals between wide field of view channel (PMT3) and the spectrometer channels 
(PMT1 and PMT2). Since the polarization of the transmitted laser pulse is rotated by 
90° between laser pulses, each channel alternatively receives perpendicular and parallel 
components. The received backscatter signal is prefiltered with an interference filter and a 
low resolution etalon pair before being directed into a beam splitter. The signal detected 
with PMT1 contains the information about the total aerosol and molecular backscatter 
signal. The signal directed through the iodine cell and detected by PMT2 is a combination 
of a amount of aerosol backscatter signal which passes through the absorption cell and the 
wings of the molecular backscatter signal. 
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After background filtering, the signal is divided into two channels by a beam 
splitter. The signal detected with the PMTl (combined aerosol+molecular channel) 
contains the information about the total aerosol and molecular backscatter signal. 
The signal directed through the 43 cm long iodine cell and detected by the PMT2 
(molecular channel) is a combination of the small amount of aerosol backscatter signal, 
which passes through the absorption cell, and the wings of the molecular backscatter 
signal (see Figure 7). The iodine absorption cell is described in more detail in Chapter 
4. 

A beamsplitter with a 30/70 (trans./refl.) splitting ratio is used to divide the 
signal between PMTl and PMT2. The beam splitting ratio is chosen from the com- 
mercially available beamsplitters so that a good photon counting statistics with both 
channels is achieved when the system is tuned to the iodine absorption peak. Some 
beamsplitter types have poor transmission and reflection efficiencies and therefore 
part of the signal is lost in the beamsplitting process. Because good system transmis- 
sion is important, a beam splitter with a high transmission and reflection efficiencies 
was chosen. 

The beam splitters are highly polarization sensitive elements. When linearly po- 
larized laser beam interacts with atmospheric aerosols and molecules, the backscatter 
signal consist of two perpendicularly polarized components. Therefore, the polariza- 
tion sensitivity of the beamsplitter transmission and reflection has to be taken into 
account. The polarization sensitivity of the beamsplitter performance is noted as a 
problem in a system described by Krueger et a/. 17 . The advantage of the HSRL is 
that the received signal is polarization filtered and therefore the light coming from 
the atmosphere to the beam splitter is always incident with the same polarization. 
On the contrary, the light from the current calibration light source is polarized or- 
thogonally to the lidar signal and therefore the polarization sensitivity of the beam 
splitter creates a problem. This problem is solved by generating a circularly polarized 
light for the beamsplitter with a ^-waveplate. 
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Figure 7. The spectral response of the HSRL receiver. The insert figure at the upper right 
corner shows a simplified setup of the spectrometer channels. The light that passes the 
field stop contains information about the Doppler- broadened molecular spectrum and the 
aerosol backscatter spectrum. The signal is divided into two parts with the beam splitter. 
The signal, which is transmitted through the beam splitter is detected with PMT1 and it 
contains information from the total backscatter signal. The signal, that reflects from the 
beam splitter is directed through an iodine absorption cell. In t his process, the aerosol 
backscatter signal is attenuated in the iodine absorption line. Also the central parts of 
the molecular spectrum are filtered out and therefore the signal detected with PMT2 is a 
combination of the wings of the molecular spectrum and that part of the aerosol backscatter 
signal that passes through the cell. 
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3.3.2 Data system 

In order to increase the dynamical range of the HSRL, new photon counting electronics 
were designed (Figure 8). The photon counting data system has three counters and 
each counter has 2 buffers, one for each polarization. Counting rates have been tested 
up to 350 MHz, but counters are designed to operate up to 1 GHz. The 100 ns bin 
width enables a 15 m range resolution. The 4 kHz laser repetition rate enables the 
use of 2333 range bins, which can be used to probe the atmosphere up to 35 km. The 
limiting factor of the data system is the photomultipliers (EMI 9863B/100), which 
allow counting only up to 15 MHz without significant pile-up effects. The PMT’s are 
pile-up and afterpulse tested and both corrections are included to the data processing. 
The corrections are discussed later in this chapter. 

Computer control of the system is made by using Intel 80960CA embedded pro- 
cessor. The user interface is realized by using a Sun Sparc Station II. In addition to 
the data taking and the communication with the Sun, the system processor is used 
to control HSRL operation including pressure controlling of the etalons, timing of 
the Pockels cell voltage, controlling of the seedlaser temperature, controlling of the 
system shutters and relays, and controlling of the WFOV aperture. 

Most of the HSRL-operations are computer controllable from a menu interface. 
System performance can be followed from real time displays. The real time Range 
Time Indicator (RTI) displays for the raw and/or inverted data are used to display the 
incoming data so that altitude, variability, thickness and depolarization of the cloud 
layer can be observed when the clouds move over the lidar site. Similar displays for 
the cloud optical properties (optical depth, aerosol backscatter cross section, phase 
function) can be shown. 
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Figure 8. The HSRL data acquisition system. 

The new counter design has improved the background light measurement accuracy 
of the HSRL. Instead of interrupting the data taking every 3.5-7 min for background 
and calibration checks (as was made with the system described by Grund 10 ), the new 
HSRL is capable of measuring the background light simultaneously with the data. 
The maximum resolution of the data system is increased from 256 to 4096 range bins 
per channel. In the HSRL, the amount of background light in the measurement is 
calculated from the measured profile by using the upper altitudes of the lidar profile. 
The signal together with the background light is saved and therefore different kind of 
background corrections can be studied. Currently, the measured HSRL profiles are 
background corrected by assuming the signal above ~33 km to be background. By 
averaging the upper range bins together, the average amount of background in the 
measurement is calculated. The disadvantage of the current background correction 
method is that some of the signal can still be originated from Rayleigh scattering 
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leading into an overestimate of the amount of background light for altitudes 33 km or 
higher. A development of a background correction method that removes the Rayleigh 
scattering contribution from the background is currently in progress. 

Photomultiplier afterpulsing can be seen as multiple output counts for a single 
input photon 30 . The afterpulsing could be caused by positive electron production 
in the dynode chain. The afterpulse pulses are time delayed compared to the signal 
from incident photon and the time delay is proportional to the path length for electron 
propagation. The signal originated from a short light pulse and measured with the 
HSRL can be presented in following form 

yf-f df 

< S(t) >= jf S(t')*(( - t')M, (14) 


where 

< S(t) > = signal measured with the HSRL 

S(t’) = signal incident on the receiver 

$ = signal responce function, that contains the laser pulse width, 

datasystem bin width and afterpulsing of the photomultiplier 

The afterpulse probability distributions were measured by illuminating the photo- 
multiplier photocathode with a short duration laser light pulse. Pulse counting rates 
less than one count per 100 ns long data bin per pulse were used. A large number 
of laser pulses were averaged (-1 hour average). The probability distributions were 
calculated from the background corrected data. The amount of background light was 
calculated from the end of data record and the average value was subtracted from 
data. The measured afterpulse probability distributions of two photomultiplier tubes 
are presented in Figure 9. The afterpulse probability distribution was found to be 
different for each tube and therefore a separate correction for each tube had to be 
implemented to the signal analysis program. With a known afterp uls e probability 
distribution function, the afterpulse contribution can be removed from the measured 
signal by using deconvolution. The effect of afterpulse correction to the measured 
signal is presented in Figure 10. 
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Figure 9. The temporal distributions of two EMI 9863B/100 photomultiplier tubes. The 
temporal distribution of the measured photons is a convolution of the laser pulse width, 
bin width of the data system, and the afterpulsing of the photomultiplier. The laser pulse 
is located in the first bins and the afterpulse can be seen as a decaying tail. The afterpulse 
distribution of the PMT1 shows better afterpulse behavior than the PMT2. 
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Figure 10. The effect of afterpulse correction tn 

thick water cloud (at 5.6 km) with an ire r i measured aerosol signal in the case of 
increases .he aerosol signal “ 51 afterpidsing 

Of the clouds is underestimated. The effect of afte ’ 7‘ the geometrica l thickness 

and thus a slight overestimate of the f * t0 ^“ ol -ular signal is similar, 
correction is not performed. P °* ^ doud 1S made > if the afterpulse 


and tkeyZ^T “° ^ **« «« «* distinguishable 

. . 118 Causes a a0nlmear response. In order to be able to separate 

between counts, a minimum separation time bet^en two puises is retired “ 

e distinguished. This time is called resolving time (or dead time) Since the 

une mterval distribution o, photons that reach the photomultiplier follow a PoLon 
distribution, the Drobabilitv p *u. , Poisson 

certain time to, V • ° * ? "**** ^ another P ulse ^ a 
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Po = 1 - e~ Rr 


( 15 ) 
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where 

r = resolving time 

R = average count rate 

In order to obtain the resolving time (dead time) of the photomultipliers, the 
photomultiplier responses to different count rates were tested. First, light level is set 
to level which is high enough to generate a pile-up. Then the photomultiplier response 
to different light levels is measured by changing the light intensity with neutral density 
filters and recording the signal change. In order to be able to define the resolving 
time of the photomultiplier, the photon counting was modeled. The pulses coming 
from the photomultiplier were treated as Poisson distributed signals. The resolving 
time that produced the best fit between simulated and measured results was used 
for the pile-up correction. All photomultipliers were found to have ~ 13 ns resolving 
time. The pile-up correction has to be performed before the afterpulse correction or 
any other corrections, because the pile-up effect is nonlinear. The pile-up is affected 
by the signal strength, signal background and afterpulsing and therefore, data has 
to contain the information from all these factors when the correction for pile-up is 
performed. 

In order to make the pile-up correction for data, the following equation was used 

N = N 0 e~ RT (16) 


dt 

This equation is solved iteratively for N 0 , so that 

1. N 0 ( 1) = N 

2. New N a (i + 1) = Ne 2i >r , where dt = bin length [s] 


where 

N = measured counts 

N 0 = actual counts 

r — resolving time 

R = average count rate = 
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3. Step 2 is repeated until N a converges 

The effect of pile-up correction to the measured signal is presented in Figure 11. 



Figure 11. The effect of pile-up correction to the measured aerosol signal in the case of 
thick water cloud (at 5.6 km) with an ice crystal precipitation at 5.1 km. The calibration 
fiber signals are shown as spikes at altitudes ~ 0 km. The maximum number of photons 
measured for the water cloud is 0.68 photons/100 ns bin and after pile-up correction 0.75 
photons/ 100ns bin are observed. Even with this low count rates the effect of pile-up is ~ 
10% and therefore, the pile-up correction is very important when strong returns from clouds 
are studied. 
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4 Iodine absorption filter 

In order to overcome the limited capability of the high resolution etalon to separate 
between aerosol and molecular scattering and to increase the system stability and 
reliability, an iodine absorption filter was constructed. For the first HSRL measure- 
ments a 43 cm long cell was made. The cell was made from glass tubing with an 
attached side arm. Optical quality end windows with anti-reflection coatings were 
epoxed to the ends. The cell with iodine crystals in a side arm was evacuated and 
kept at 27 °C. Transferring of the iodine from the absorption cell into the vacuum 
pump was prevented by evacuating the cell through a cold trap and cooling the side 
arm with liquid nitrogen. Although the iodine cell can be operated at room temper- 
ature, the operating temperature of the cell has to be controlled, because the vapor 
pressure of iodine is very temperature sensitive 32 . In the HSRL, the cell temperature 
is maintained with ±0.1 °C accuracy by operating the cell in temperature controlled 
environment. 

The iodine spectrum is measured by scanning the laser wavelength by changing 
the temperature of the seedlaser under computer control. A small amount of laser 
light is directed into a 100 m long fiber optic delay (Fiber 1 in Figure 4) and sent 
to the receiver to create a calibration light source. The temperature-wavelength 
dependence of the scan was determined by using the free spectral range of the high 
resolution etalon as a reference. This could be made, because the free spectral range 
of the etalon can be calculated when the length of the etalon spacers is known and 
the spacing of two (or more) etalon transmission peaks in temperature units can be 
measured. The calibration was made by simultaneously measuring the transmission 
spectrum of the high resolution etalon and iodine absorption filter. The simultaneous 
measurement of the high resolution etalon and iodine absorption filter transmissions 
was made by measuring the signal reflected from the high resolution etalon (Figure 
12) and the signal transmitted through the absorption cell. The pressure in etalon was 
held constant while the laser wavelength was scanned. The spectrum was normalized 
by measuring the cell transmittance without the iodine cell. 

A part of the measured iodine spectrum is presented in Figure 13. The measured 
spectrum was compared with a published spectrum 18 and an ~0.01 pm wavelength 
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agreement in relative line positions was observed. The linearity of the temperature 
scan was confirmed from the free spectral range information of high resolution etalon 
by performing the scan over more than one free spectral range. Single mode operation 
between two seedlaser mode hops can be maintained over 20 GHz range (at 1064 nm) 
and within this range two high resolution etalon free spectral ranges can be covered. 
During a mode hop the laser frequency jumps back about 10 GHz. 



Light from the 
calibration fibers 


Figure 12. The HSRL receiver used for iodine spectrum calibrations. The same receiver 
setup was used for the first HSRL measurements with the iodine absorption filter. For data 
taking the transmission of the high resolution etalon was tuned out from the peak and the 
etalon was used as a reflector. When the beamsplitter and the mirror 4 axe removed from 
the system, the system returns back to the old HSRL receiver. 
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Figure 13. Transmission of a 4 cm and 43 cm iodine cells as a function of wavelength shift. 
The identification fine numbers are from Gerstenkorn and Luc 18 . 

For initial HSRL measurements the line 1109 (peak wavelength 532.26 nm) 18 , 
which is well isolated from the neighboring lines, was chosen. The full width half 
maximum width of the line is ~ 1.8 pm and the peak tr ansmiss ion is ~ 0.08%. The 
hyperfine structure of the peak 1109 defines the asymmetric shape of the absorption 
peak 33 . In fact, the line 1109 is a combination of two rotational vibrational transitions 
with different hyperfine structures. 

The iodine absorption cell provides a robust filter for the HSRL, because it is not 
dependent on the mechanical alignment of the filter or the angular dependence of the 
incoming light. Another advantage is the stability of the absorption characteristics. 
This provides a stable long term operation. The strength of observed absorption line is 
dependent on the line strength, and the length, temperature, and pressure of the cell. 
By controlling the operating environment and with a nearly leak proof system, the 
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current iodine cell is operated for several months without any maintenance. During 
this time, a small change in absorption strength and line width were observed due 
to a small leak that was caused by the iodine penetrating through a hose. Also the 
iodine was found to condense into the walls of the cell, but even during a long period 
of time, the amount of condensation has been small and ~ 10% extra absorption 
is observed. The condensation can be prevented by operating the tip of the side 
arm couple degrees below the cell temperature. The problems with reactive iodine 
penetrating through the hoses can be prevented by using a sealed all-glass cell. In 
a short term operation, the stability of the absorption characteristics has proven to 
be so good that a system calibration scans from different days can be used for the 
calculations of the system calibration coefficients. This requires, that the alignment 
of the receiver optics is stable. 

An absorption filter offers a high rejection against aerosol scattering and therefore 
it makes the separation between aerosol and molecular scattering easier. Also, a wide 
dynamic range in rejection against aerosol scattering is achieved by simply changing 
the vapor pressure or the length of the cell. Comparison between high resolution 
etalon and iodine absorption filter performance is presented in Figure 14. A 2:1 
separation between molecular and aerosol scattering by the etalon (Figure 14. b) is 
measured compared to a 1000:1 separation in the iodine cell when operated at 27 °C 
(Figure 14.a). The molecular transmission in Figure 14.a and Figure 14.b is calculated 
by using the Doppler- broadened molecular spectrum at -65 °C. This temperature is 
close to the lowest temperature measured at the tropopause and this gives the smallest 
transmission through the iodine absorption cell. The molecular transmission of the 
high resolution etalon and the iodine absorption filter are similar (Figure 14.c). Due 
to wide absorption line width, the molecular transmission of the iodine filter is more 
dependent on the air temperature than the etalon. The temperature dependence of 
the cell transmission is modeled by using the table values of iodine vapor pressure^ 
(Figure 14.d). Calculations show, that by changing the cell temperature from 27 °C 
to 0 °C, the online transmission can be tuned from 0.08% to 60%. 



37 



Wavelength Shift (pm) 



- 3 - 2-1 0 1 2 3 

Wavelength Shift (pm) 


c 

o 

To 


CO 

E 


CO 

c 

CO 


K 



Figure 14. (a) Transmission of 43 cm cell (solid line) together with the molecular transmis- 
sion (dashed line) at -65 °C air temperature as a function of wavelength shift. Dot-dashed 
line shows the calculated molecular spectrum at -65 °C. (b) Etalon tr ansmissi on (solid line) 
and calculated molecular transmission (dashed line) as a function of wavelength shift. Dot- 
dashed line shows the calculated molecular spectrum at -65 °C. (c) Comparison of molecular 
transmission of high resolution etalon and iodine cell as a function of air temperature, (d) 
Iodine cell aerosol and molecular transmission as a function of cell temperature. 
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5 Calibration and tuning 

5.1 Calibration theory 

The HSRL measures optical properties of aerosols by using the backscatter from 
atmospheric molecules as a calibration target. In the receiver the signal is separated 
into two separate signals: one proportional to total aerosol and molecular scattering 
and the other containing the molecular backscattering together with a small aerosol 
cross talk component. The measured signals are 

S a+m (R) = 7 / [N a (R) + N m (R)} ( 17 ) 

Sm(R) = V [CamN a (R) + C mm (R)N m (R)\, ( 18 ) 

where 

S a +m = signal measured with the combined channel (PMT1 in Figure 6) 

S m = signal measured with molecular channel (PMT 2 in Figure 6) 

N a ,N m — total number of aerosol and molecular backscatter photons 

incident on the receiver field of view 
Cam = aerosol transmission of the molecular channel relative to 
the combined channel 

Cmm = molecular transmission of the molecular channel relative to 
the combined channel 

77 = system efficiency factor that includes the optical transmission 

of the combined channel and its photomultiplier quantum efficiency 
These two equations can be solved to present the separated aerosol and molecular 
backscatter signals. 


N m (R) = 


(R) Cam ‘S'o+m ( R) 

V(C mm {R) ~ Cam) 


(19) 


N (R) = V N m(R) 


( 20 ) 


The calibration coefficients C am and C mm are obtained from a system calibration 
scan. For calibration the system input aperture is uniformly illuminated with a diffuse 
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light. The receiver spectral transmission function is measured by scanning the laser 
wavelength over an 11 pm wavelength range around the selected iodine absorption 
peak and recording the signals (originated from calibration fiber 1 and calibration fiber 
2) with both spectrometer channels (PMT1 and PMT2 in Figure 6). A calibration 
scan is performed before and after each dataset. When system is running for a long 
period of time (time ~ 3 hours) the system operation is interrupted and a calibration 
scan is performed. An example from a calibration scan is presented in Figure 15. 
In addition to the information on the system spectral transmission, the calibration 
signals contain information on the beamsplitting ratio between channels. Since the 
signal measured through the iodine cell is flat at the top of the iodine absorption 
peak, the determination of the wavelength of the absorption maximum is based on 
the signal from the reference iodine cell (4 cm cell in Figure 5, calibration fiber 2) 
and measured with the PMT1. Otherwise, the signals from the first calibration fiber 
are used for the calibration coefficient calculations. 



Wavelength Shift (pm) 


Figure 15. An HSRL calibration scan. The calibration fiber 1 signal that is detected with 
PMT2 shows the iodine absorption spectrum of the 43 cm long iodine cell. The calibration 
fiber 2 signal detected with PMT1 presents the absorption spectrum of the 4 cm long 
reference cell. The signal from calibration fiber 1 and detected with PMT1 is used as a 
reference. 
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Since the Doppler-broadening of the aerosol backscatter is negligible, the spectral 
distribution of the aerosol backscatter can be assumed to be similar to the spectral 
distribution of the transmitter laser. The measured calibration signals can be pre- 
sented as a convolution between laser spectral distribution and spectral bandpass of 
each channel. Therefore, the fraction of the total aerosol backscatter detected by the 
molecular channel (C om ) can be directly obtained from the calibration signals. 


r _ S m (at peak) 

° m S a+m (at peak) ’ ( 21 ) 

where 

S m = calibration fiber 1 signal detected with the PMT2 at the 
iodine absorption peak 

S am = calibration fiber 1 signal detected with the PMT1 at the 
iodine absorption peak. 

The fraction of the total molecular backscatter measured by the molecular channel 
(Cmm) is calculated by convoluting the measured filter function with the calculated 
molecular spectrum. 


, 5'm(-^n)M a pect(A n )AA 

Sn=l (\ ) Mipect ( \ ) A A 


( 22 ) 


where 

S m = calibration fiber 1 signal detected with the PMT2 

(filter function for molecular channel convoluted with the 
laser spectrum) 

Sa+m = calibration fiber 1 signal detected with the PMT1 

(filter function for aerosol + molecular channel convoluted with the 
laser spectrum) 

Aspect — calculated molecular spectrum 

N = number of points in calibration scan 

A = wavelength 
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AA = the wavelength difference between two points in the calibration scan 

The divisor on the Eq. 22 is presented as a convolution aerosol and molecular 
channel. Therefore, the divisor presents the amount of molecular spectrum seen with 
the combined aerosol and molecular channel. The dividend of the Eq. 22 describes the 
molecular signal detected through the iodine absorption cell. The molecular spectrum 
model used in the calculation is presented in a paper by Yip and Nelking 34 and it 
includes the effects of Brillouin scattering as a function of temperature and pressure. 

The accuracy of the calibration coefficients is mainly limited by the photon count- 
ing statistics. Because the signal transmitted through the absorption peak is small, 
the error due to photon counting statistics dominates the error in the determination 
°f Cam- Therefore, the accuracy of the C am is improved by increasing the photon 
counting statistics at the absorption peak. Three different ways to increase the pho- 
ton counting statistics can be considered. First, the signal at the absorption peak can 
be increased by scaling the light with neutral density filters while scanning. Second, 
the amount of aerosol backscatter signal can be further decreased into a point where 
the effects of the photon counting statistics are negligible. Third, longer averaging 
time can be used. 

The disadvantage of using neutral density filters is that the filters have to be well 
calibrated and the change in the value of neutral density filter has to recorded into the 
data so that the signal can be reconstructed back to the absorption spectrum. The 
disadvantage of the longer absorption cell is that the increased cell length will further 
decrease the amount of transmitted molecular signal. Also the spectral purity of the 
laser limits the observable absorption strength. In order to be able to obtain a good 
photon counting statistics for the signal of the whole absorption peak, a long averaging 
time is required and therefore the total calibration time would be unreasonable long 
(~ It) and during this time the laser has time to drift. The drift in the laser output 
wavelength during the scan effects the width of the measured absorption spectrum. 

The current HSRL uses a calibration procedure, where the absorption spectrum 
is first measured by scanning the laser wavelength so that ~ 1% photon counting 
accuracy is achieved for the spectrum around the absorption peak. In order to obtain 
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a high photon counting statistics in short period of time, the light from the calibration 
fibers is optimized so that maximum number of photons is detected with small pile-up 
effects at the detectors. During the scan the location of the peak absorption maximum 
is detected from the signal through the 4 cm long reference iodine absorption cell. 
After completing the scan, the seedlaser temperature is set back to the maximum 
and by using a tuning program (described in more detail in Chapter 5.2) the laser 
wavelength is kept at the absorption peak until better than 3% photon counting 
statistics is obtained. With this procedure the effects due to a shift in the laser 
output wavelength to the width of the absorption spectrum can be minimized and 
the photon counting errors in the determination of the C am can be reduced from 
about 20% to 3% within ~ 10 min averaging time. 

The atmosphere provides the best reference when the accuracy of the HSR1 cal- 
ibrations is studied. Figure 16 presents an HSRL calibration which is performed 
simultaneously with data taking. Two different cases are studied. First, a calibration 
from a thick water cloud is shown. Second, a calibration from clear air is presented. 
In order to detect the possible range dependence of the calibration, lidar returns 
from different altitudes are studied. The comparison between calibrations from at- 
mosphere and from the calibration light source also recovers possible misalignments 
of the system. 

The system calibration signal from the iodine absorption spectrum presents a cal- 
ibration from a pure aerosol target. The agreement between system calibration and 
atmospheric calibration from a thick water cloud can be seen from Figure 16. Both 
signals are defined from the ratio of the signal detected through the iodine cell to 
the signal detected with the combined aerosol+molecular channel. The background 
corrected, energy normalized signals are used. The data is averaged over a 90 m 
range. An expected calibration curve from a pure molecular target can be calcu- 
lated by convolving the measured iodine absorption spectrum with the calculated 
molecular spectrum. The calculated molecular calibration together with a measured 
atmospheric calibration from different altitudes are presented in Figure 16.b-d. The 
measured absorption spectrum is presented as a reference. For the calculated molecu- 
lar calibration, the atmospheric temperature, and therefore the width of the Doppler- 
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broadened molecular spectrum, is calculated by using the temperature values obtained 
by a radiosonde measurement. The signals from higher altitudes are disturbed by the 
low photon counting statistics, but otherwise a good agreement between system cal- 
ibration and atmospheric calibration is obtained and no range dependence in the 
system calibration is observed. The range dependence of the atmospheric calibration 
would show up as a noticeable deviation from the system calibration. 
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Wavelength Shift (pm) 

Figure 16. A HSRL calibration scan together with a simultaneous calibration from the 
atmosphere. Figure (a) shows a calibration from a thick water cloud (thin dashed line) 
together with a system calibration scan (thick solid line). In figures (c)-(d), the dashed line 
shows a clear air calibration at 3175 m (b), 5510 m (c), and 7550 m (d). The temperatures 
at these altitudes were -11 °C (b), -32 °C (c), and -45 °C (d), respectively. The long dashed 
line presents the expected molecular return. The measured calibration scan is presented as 
a reference (solid line). 
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5.2 The laser wavelength locking into the iodine absorption 
peak 

The system calibration is sensitive to the drifts between the transmitter wavelength 
and the receiver bandpass transmission maximum. The measurements show that 
the seedlaser drifts at ~100 MHz/h rate. In order to achieve a stable long term 
operation without frequent calibrations, the wavelength of the transmitter laser is 
locked to the iodine absorption peak. Wavelength locking to the Doppler-broadened 
iodine absorption line was used because it requires much less power than a locking 
into a hyperfine structure of the iodine absorption line with Doppler- free technique 33 . 
Compared to the wavelength locking with a high resolution etalon, the advantage of 
the locking to an iodine absorption peak is that, iodine peak provides an absolute 
frequency reference. Another technique to lock the laser output wavelength to the 
Doppler-broadened iodine absorption line was reported by Arie and Byer 35 . They use 
Fourier transformation spectroscopy to lock the laser to the center of the Doppler- 
broadened peak. This method does not require any dither of the laser frequency, but 
it is more complex. 

In the HSRL, the absorption spectrum of the 4 cm long iodine absorption cell is 
used to provide information about the absorption peak maximum. The absorption 
peak of the 43 long iodine absorption cell cannot be used as a reference for the 
wavelength locking because the saturation of the absorption at the peak causes the 
flat shape of the peak, and because the signal at the peak is small due to the strong 
absorption. Therefore it does not provide good photon counting statistics for the 
locking. The length of the reference cell is chosen so that the absorption is — 50%. 
The cell transmission has to be high enough to provide a good photon counting 
statistics within a short averaging time. 

The locking of the laser wavelength to the iodine absorption peak is performed 
by using an automatic controlling program that works as follows. First, the location 
of the absorption maximum is detected during a calibration scan. After completing 
the scan, the program automatically sets the seedlaser temperature to the observed 
peak. In order to keep the laser wavelength locked to the maximum absorption 
wavelength, the seedlaser temperature is dithered around the optimum temperature 
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and information about the ratio of signal from the second calibration fiber to the signal 
from the first calibration fiber detected with PMT 1 (see Figure 17) is gathered. The 
seedlaser temperature is kept at temperature that produces the minimum ratio. The 
basic idea of the tuning program is presented in the following. 

The tuning procedure has four steps. First, information about the ratio be- 
tween calibration fibers is gathered for the seedlaser temperature (T(peak)) that 
was detected to produce the minimum ratio. Then, a temperature T=T(peak)+dT 
is applied and the change in the ratio is observed. After gathering enough statistics 
(~ 30 s), the seedlaser temperature is returned back to the temperature T(peak) and 
a new value for the ratio at this temperature is measured. After this, the optimum 
temperature is detected by finding the temperature that produces the minimum ratio. 
If the temperature T=T(peak)+dT produced a smaller ratio, then that temperature 
becomes to the new optimum temperature T(peak). If the ratio for temperature 
T=T(peak)+dT was not better, the dithering to temperature T=T(peak)-dT is per- 
formed and the procedure is repeated. The wavelength dither corresponds to temper- 
ature change of 0.009 °C (0.052 pm). A time history of the dither temperatures for a 
9 hour run shows that, the seedlaser temperajture is dithered between 3 temperatures 
under typical operating conditions. 
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Wavelength Shift (pm) 


Figure 17. Signal for the frequency locking of the laser. An expanded view from the peak 
shows that as the temperature is changed in either direction from the detected peak, a 
change in the ratio is noticed. 

Because the absorption peak of the 43 cm long absorption cell is flat around the 
maximum absorption wavelength, the error due to tuning to the amount of aerosol de- 
tected in the molecular channel (C om ) is comparable to the photon counting statistics. 
On the other hand, the error due to dithering to the amount of detected molecular 
signal in the molecular channel ( C mm ) is a combination of photon counting statistics 
and the error between the convoluted signals at different dithering wavelengths. The 
errors in C mm due to tuning as a function of atmospheric temperature are presented 
in Figure 18. Figure 18 shows the error in the determination of C mm , when the laser 
wavelength is tuned off by ±0.052 pm, but the calculation of C mm is made for the 
peak wavelength. Also the error due to dithering is shown. The asymmetry of the 
absorption spectrum makes the errors due to the tuning asymmetric. The total effect 
of the tuning procedure to the measured profiles has to be calculated as a weighted 
average of the errors at different dithering wavelengths, because the tuning program 
is realized so that the laser spends 2/3 of the time at the wavelength that produces 
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the minimum ratio and 1/3 of the time doing the dithering. Therefore, the total error 
due to dithering is better than 0.1%, when measurement period is long compared to 
the dither time. In principle, the error due to wavelength locking can be eliminated 
by inverting the data by using different calibration coefficients for different dithering 
wavelengths. This has not been accomplished yet. 
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Figure 18. The errors in C mm due to the tuning. Error when the seedlaser temperature is 
detuned from the optimum temperature by ± one step (0.009 °C) and when the seedlaser 
temperature is dithered, but the inversion is performed by using the observed peak value. 
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6 Measurements 

Starting from July 1993, the iodine absorption filter based HSRL has been operated 
at the University of Wisconsin-Madison campus. During this time the HSRL has been 
routinely operated and the stability and reliability of the system have been tested. 
As a result, a dataset consisting of ~ 30 different cirrus cloud cases has been obtained 
with a simultaneous NOAA-11 and/or NOAA-12 satellite overpasses. As an example 
from the data collected with the HSRL a dataset from November 11, 1993 is shown. 
This dataset contains a meteorologically interesting case: a cirrus cloud, supercooled 
water cloud, and ice crystal precipitation together with a strong low level aerosol 
structure. Figures 19 - 25 present Range Time Indicator (RTI) pictures from the 
data. Both raw and inverted data are shown along with depolarization and optical 
depth. The pictures are generated from the background corrected, energy normalized, 
and range square corrected data. 

For the RTFs of the inverted data, the aerosol and molecular signals are separated 
by using Equations 19 and 20. The optical depth is obtained from the ratio of the 
inverted molecular profile to the return predicted for the pure molecular scattering 
(Eq. 10). The color scale shows the signal strength and the white areas are regions 
where the backscatter signal is larger than the maximum color scale value. The black 
areas indicate that the signal is smaller than the smallest color scale value. 

The Figure 19 shows the raw lidar return detected with the combined aerosol and 
molecular channel. This profile is similar to the profile obtained with a conventional 
single channel lidar: the signal from small amounts of aerosol scatterers is dominated 
by the scattering from molecules, and therefore all aerosol structures are not clearly 
visible. The ability of the HSRL to separate aerosol and molecular scattering can be 
seen from the RTI picture of the inverted aerosol signal Figure 20. After inversion, 
the aerosol structures are more visible and they do not have the decrease with altitude 
caused by the atmospheric density profile. The ability of the iodine absorption filter to 
reject aerosol scattering is visible from the RTI of the raw molecular signal (see Figure 
21): only a small aerosol cross-talk for the densest parts of the clouds is observed and 
this is easily removed by the inversion, as can be seen from Figure 22. The phase 
of the water at different layers can be seen from the depolarization RTFs. Figure 
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23 shows the raw depolarization observed with the combined aerosol and molecular 
channel. The inverted aerosol depolarization is shown in the Figure 24. From these 
pictures, a cirrus cloud at ~ 8 km (depolarization ratio ~ 40%) and a supercooled 
water cloud at 5 km ( depolarization ratio ~ 1 %) with ice crystal precipitation can be 
easily separated. For low level aerosols (0-3.7 km), a two layer polarization structure 
is seen. The small increase in water cloud depolarization as a function of cloud height 
is an indicator of multiple scattering. The low molecular depolarization is presented 
in Figure 25. The low depolarization ratio values with small signal to noise ratio show 
up in the picture as noise. 

The optical depths on the different parts of the data set can be seen from the 
Figure 26. The optical depth of the cirrus cloud between 7 and 10 km is ~ 0.4. The 
water cloud at ~ 5 km has an optical depth of 2.5-3. The extinction through the 
ice crystal precipitation below the water cloud and the extinction through the water 
cloud can be seen as a change in the color scale as a function of altitude. The optical 
depth of the ice crystal precipitation is ~ 0.1. 

A more detailed analysis of the dataset is presented in the following sections. First, 
the depolarization measurements are discussed in Chapter 7.1. The effects of mul- 
tiple scattering to the depolarization measurements are shown. The measurements 
of the cloud particle sizes are not included to this study. The depolarization data 
from August to November 1993 is analyzed and a summary from the observed depo- 
larizations as a function of atmospheric temperature is given. Second, an example 
from a measurement of scattering ratio, aerosol backscatter cross section, and optical 
depth is given together with error estimates for the optical depths (see Chapter 7.2). 
The temperature dependence of the Doppler- broadened molecular spectrum enables 
the measurements of the atmospheric temperature by the HSRL. Preliminary results 
from a temperature measurement are presented in Chapter 7.4. 
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Figure 19. The raw lidar return presenting the combined aerosol and molecular channel 
return. A water cloud layer with an ice crystal precipitation are seen at 5 km. Above the 
water cloud, a cirrus cloud can be seen. The low level aerosol structure between 0 and 3.7 
km is hardly visible because it is damped by the molecular backscatter signal. 
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Figure 20. The inverted aerosol profile. After separating the aerosol and molecular ba cks eat- 
ter returns, the layers where the aerosol backscatter signal is small compared to the molec- 
ular signal are clearly visible. The largest difference is seen for the low level aerosol layer 
between 0 and 3.7 km. > • 




Figure 21. The raw molecular return. A small aerosol cross talk signal is visible, for the 
densest parts of the water cloud at ~ 5.5 km. The dark areas indicate that very little, pr no 
return through parts of the water cloud is observed. 
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Figure 22. The inverted molecular profile. After the inversion, the cross talk ;t)bat yas 
visible in Figure 21 cannot be seen and the inverted molecular profile therefore- presents 
the atmospheric extinction at various points of the dataset. The inaccuracy of the overlap 
correction can be seen as a darker line at ~ 1 km. 







55 



Figure 23. The raw aerosol depolarization combined with the molecular depolarization. The 
picture shows the depolarization ratio that is seen with a lidar that cafinot separate the 
aerosol and molecular backscatter signals. The depolarizations for altitudes with low aerosol 
content are dominated by the molecular depolarization. The parts of thd^irrus clcKid and 
parts of the ice crystal precipitation between 4 and 5.5 km show depolajii&$^ 
are ~ 10%, and those layers could be expected to contain mixture of ice and ' water. Some 
parts of aerosol layer between 2 and 3.7 km show depolarization of ~3.5 %. The water 
cloud at 5.5 km has ~1% depolarization. The increase in the water cloud depolarization as 
a function of altitude is due to the multiple scattering. 



.iSfrri*- 





56 






o -t- 


; ' * 


?r —c-u- 

^ ;;>. ■■ •..•■:■ , : v 


]^s--i£k 



>*vV 1 

... j'f^yU^ 

Figure 24. The inverted aerosol depolarization. After inversion, the y^s cloud 
ization of ~ 40% indicates pure ice depolarization. The ice crystal predpitatibn 
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from the water cloud show similar depolarizations values. The low levei^aerosol structurG 
shows a two layer polarization structure. A ~1% depolarization ratio foiMtiie layer between 
0.5-2 km is observed indicating nearly spherical particles. The depolarization of the ; layer 
between 2 and 3.7 km shows a ~ 5% depolarization. * 
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Figure 25. The inverted molecular depolarization. A less than 1 % 
is observed. The increase in the depolarization variations at the 
the low signal to noise ratio. 
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Figure 26. The optical depth. The optical depth above 3 km is shown. The optical depth 
of the cirrus cloud between 7 and 10 km is ~ 0.4. The water cloud at 5 km has a" optical 
depth of 2.5-3. The extinction through the ice crystal precipitation below the water cloud 
and the extinction through the water cloud can be seen as a change in the color scale. The 
optical depth of the ice crystal precipitation is ~ 0.1. 




6.1 Depolarization ratio 
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Lidar measurements of atmospheric depolarization can be used to distinguish between 
liquid and solid phases of water in the atmosphere 36-38 . The quantity to describe the 
degree of polarization is the linear depolarization ratio S 


6 = 


’ ts ■ . v 


. : ( 23 ) 


where /j_ and /y are the measured perpendicular and parallel back'^eatter intensities 
in respect to the transmitter polarization axis. The backscatter signal of a linearly 
polarized laser beam from spherical particles is totally linearly polarized (<5 = ’6). /The 
particles can be assumed to be spherical in case of wet haze, fog, cloud droplets, and 
small raindrops. The depolarization of the pure molecular atmosphereris jipnzero, 
because of the anisotropy of the air, and therefore a 0.37-0.4 % depolarization for 
the Cabannes line is expected 39,40 . The depolarization of the molecular return that 
includes the Cabannes line and the rotational Raman lines is ~ 1.5% 40 . If particles 
are nonspherical (as ice crystals, snow flakes or dust particles) or if the backscatter 
signal has a multiple scattering contribution, the backscattered signal contains a cross- 
polarized component (0 < 8 < 1). A specular reflection from a oriented ice crystal 
layer provides a small depolarization, which therefore can be misinterpreted as the 
backscatter signal from a water cloud, but the off vertical pointing direction, of the 
HSRL is expected to prevent this. 

Lidar studies of atmospheric polarization have been traditionally based on a tech- 
nique, where a linearly polarized laser beam is sent to the atmosphero^^ea^dved 

.... , fk /..’I 

polarization components are separated by a polarization cube and detected by a pair 

a photomultiplier tubes, one for each polarization component. This method requires 
a precise calibration of the receiver in order to avoid the problems clue to diffwences 

■ ■ : i •. >J 4 1 , l y i v * jLt -■ * ; 

between channels in optics, and in photomultiplier sensitivity. 

In the new HSRL (see Figs. 5 and 6), the polarization measurements are made 
by using one transmitter laser and one detector for both polarization components. 
The transmitted laser beam is linearly polarized, and for polarization measurements 
the polarization of the transmitted laser beam is rotated by 90 degrees on alternative 
laser pulses by a Pockels cell. In the receiver, the signal parallel to the transmitter 
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polarization axis and the signal perpendicular to that are separated and cleaned by 
a polarization cube pair. Since the polarization of the transmitted laser. shot alter- 
nates between two subsequent laser shots, the same detector can accurately irleasure 
both polarization components and therefore no calibration of the reeaver is required. 

i.rWi 'M,i V* . i *V . ' 

The 250 //s time separation between laser pulses insures, that both depolarization 
components are measured from the same atmospheric scatterers. The depolarization 
variations as a function of field of view can be studied with the new WFOV-channel. 
The effects of specular reflection from a horizontally oriented crystal layer are mini- 
mized by tilting the receiver by 4 degrees from the zenith. With the HSRL, separate 
measurements of aerosol and molecular depolarizations can be made. 

Depolarization measurements are affected by presence of multiple scattering when 
a lidar is used to probe dense clouds. Because the amount of detected multiple 
scattering is dependent from the cloud particle size, range from the cloud, optical 
depth of the cloud, and the field of view of the system, typical systems withjl-5 mrad 
field of views have difficulties separating between water and ice forjdgnse clouds.. In 
order to separate between water and ice for optically thick clouds, fie HSRL uses a 
160 //rad field of view for the spectrometer channels. 7 ' 

‘ 1 \ f ’ ' ; ^ ,\ ' r *t * v. 

Figure 27 represents an example of the multiple scattering effect on the mea- 
sured depolarization. Figure 27. a. shows the measured energy 
corrected signals from a supercooled water cloud with ice precipithtidil' r Foif ' the 
WFOV-channels, 0.22, 0.65, 1.1, and 1.6 mrad field of views are used. The WFOV- 
channel measurements are then compared to the measurements of the 0.16 mrad field 
of view channels. The amount of multiple scattering in the signal is calculated from 
the ratio of the WFOV-signal to the signal simultaneously measured with the narrow 
field of view channel (see Figure 27. b.). Figure 27. c shows the depolarization ratios 
for different field of views and the Figure 27.d shows the ratio of the. WF0Y; depolar- 
ization to the narrow field of view depolarization. For the optically.ibln , fee.. costal 
precipitation layer, the change in depolarization as a function of field c^f.yie^jfs-hardly 
noticeable. The change from ice to water shows up as a drop in 
ratio. All field of views show a low depolarization for the water cloud base- but as 
soon as the signal penetrates deeper into the cloud, an increase in the depolarization 
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as a function of field of view can be observed. Inside 300 m the multiple scattering 
effects in the cloud will increase the depolarization of the larger field of views up to 
the level which is comparable with the observed ice crystal depolarization and there- 
fore, the separation between ice and water becomes impossible. The low values of 
depolarization observed with the narrow field of view channel and with the smallest 
WFOV field of view show that the narrow field of view effectively suppresses the 
effects of multiple scattering. Therefore, the HSRL signals from these .field of views 
can be used to separate between ice and water clouds even for optically thick clouds. 
On the other hand, information from depolarization and signal strength variations as 
a function of field of view can be used to verify multiple scattering calculations. 

Theoretically, the depolarization of the spherical dropplets should be zero, but the 
depolarization ratios observed for the water cloud at 5.5 km is ~2 %. In the cases of 
water clouds with an ice crystal precipitation, the non-zero values of depolarization 
ratio can be explained by the presence of ice. Because the water cloud is precipitating 
ice, there has to be some ice mixed with water in the cloud. The scattering from the 
cloud base is ~ 20 times larger than signal from ice crystal precipitation. Based on 
this ratio, ~ 2 % (~ 1/20 x ice depolarization of 40 %) depolarization for the cloud 
base depolarization can be expected due to the presence of ice crystals with ~ 40 % 
depolarization. 
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Depolarization 0epol(WFOV)/Depol(A+M) 

Figure 27. Effects of multiple scattering on depolarization in the backseatter, return from a 
water cloud at 5.5 km with ice crystal precipitation between 4.5 and 5.3"km. 1)ata obtained 
on November 11, 1993 02:04-02:14 UT. (a) Measured signals, (b) Ratios of me4sured WFOV 
signals to 0.16 mrad field of view signal, (c) Measured depolarization ratios. T?he cloud base 
depolarization of the smallest field of view at 5.4 km is ~ 2 %. (d) Ratios of the measured 
WFOV depolarizations to the depolarization of the 0.16 mrad field of view channel. 
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The capability of the HSRL to distinguish between aerosol and molecular scatter- 
ing allows separate polarization measurements. This is important when layers with 
a low aerosol content are studied. The effects of the molecular scattering on the 
observed depolarization can be seen by comparing the raw and inverted aerosol de- 
polarization ratios. For the cases where the amount of aerosol scatterers is small, the 
signal from molecular scattering dominates the depolarization picture (see Figure 23). 
Therefore, aerosol depolarizations similar to the molecular depolarization can be seen 
and also parts of the cirrus cloud and ice crystal precipitation show depolarization 
ratios which are close to the depolarization of the supercooled water or mixture of ice 
and water (green areas in the raw aerosol depolarization RTI). After inversion a clear 
difference in the depolarizations is seen: the depolarization of the low level aerosols is 
better defined and the depolarization of the cirrus shows that the cloud contains pure 
ice crystals (Figure 24). Therefore, a clear separation between ice and water can be 
based to the depolarization ratios calculated from the inverted aerosol profiles. The 
effect of the inversion to the depolarization ratio is also visible from the Figures 28 — 
29. 


The profile of the raw aerosol depolarization shows the depolarization of the com- 
bined aerosol+molecular channel and therefore, it shows the depolarization ratio 
which is seen with a conventional single channel lidar. The raw molecular depolar- 
ization contains the depolarization component of the aerosol cross talk together with 
the molecular depolarization. In the same figure, the separated aerosol and molecular 
signals are shown with the inverted aerosol and molecular depolarizations. Figure 
28 shows a two layer water cloud with an ice crystal precipitation. The peaks of the 
water clouds are observed at 5.2 and 5.6 km and the corresponding low depolarization 
values are 1-2 %. A small increase in the depolarization with penetration depth is 
observed. The values of raw depolarization of the ice crystal precipitation are close 
to those of water and ice mixture, but the inverted depolarization ratios show a clear 
ice depolarization. After inversion the increase in aerosol depolarization at 2-3.7 km 
is also very clear. The Figure 29 shows a water cloud layer at 5.4 km with a more 
dense ice crystal precipitation. With increased backscatter signal from ice, the value 
of raw depolarization ratio indicates clear ice and therefore the effects of molecular 
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depolarization do not show up so clearly, eventhough a change in depolarization ratio 
is observed after the inversion. ' ; 1 '' ’ 

The depolarization observed in the clear air is 0. 7-0.8 %, when the observations 
are made without the low resolution etalons. When the low resolution 'etalons are 
used, the depolarization is 0.55-0.6 %. The expected value for the molecular depo- 
larization of the Cabannes line is ~ 0.4 and ~ 1.44 for the Cabannes line 39 - 40 and 
rotational Raman lines 40 . The HSRL observations show a 1.5 % depolarization for the 
signal from the rotational Raman lines and Cabannes line. Because the system filter 
bandpass admits a small fraction of the rotational Raman lines (75 % depolarization) 
and simultaneously blocks part of the Cabannes line, the molecular depolarization 
value measured by the HSRL is larger than the expected Cabannes line depolariza- 
tion. The amount of transmitted rotational Raman signal is temperature dependent. 
A model to calculate the expected depolarization was written. This model includes 
the rotational Raman spectrum for nitrogen and oxygen, molecular spectrum includ- 
ing effects of Brillouin scattering, and spectral transmissions of different .filters. The 
polarization that correspond to different atmospheric temperatures Gan be calculated. 
The ratio of rotational Raman signal to Rayleigh signal is chosen so, that a 1 44 % de- 
polarization for the clear air is observed. These calculations show, that a 0.56 -0.62 % 
depolarization is expected for the case where no low resolution etalons were used. A 
0.402 -0.425 % depolarization is expected, when one or two etalon are used. The de- 
polarization observed with the HSRL are larger than the expected values The cause 
of the additional depolarization in the HSRL measurements is currently unknown. 
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Figure 28. Raw and inverted depolarization ratios of clear atmosphere with ice crystal 
precipitation and two water cloud layers above it (November 11,1993 01:35-01:47 UT). The 
upper set shows the raw aerosol and raw molecular profiles together with the depolarizations. 
The lower set shows the inverted profiles with the inverted depolarizations. The water clouds 
are observed at 5.2 and 5.6 km altitudes. The low water cloud depolarization values can 
be observed from the aerosol depolarization figures. The small increase in water cloud 
depolarization as a function of penetration depth is due to the multiple scattering. The 
weak signal from the ice crystal precipitation is visible 4.3-5. 1 km. Because of the low ice 
crystal content, the raw depolarization ratios of the ice crystal precipitation shows values 
that would indicate mixed phase, but after the inversion a clear ice crystal depolarization 
is visible. The clear air aerosol depolarization is less than 6%. The observed molecular 
depolarization is 0.8 %. 
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Figure 29. Raw and inverted depolarization ratios of a thin cirrus cloud and a thin water 
cloud with an ice crystal precipitation (November 11,1993 01:55-02:01 UT). The upper set 
shows the raw aerosol and molecular profiles along with the raw depolarization ratios. The 
lower set shows the inverted profiles with the inverted depolarization ratios. The water 
cloud at 5.5 km has a low depolarization value ~1.5% at the cloud base and it increases 
towards the cloud top due to the multiple scattering. The dense ice crystal precipitation is 
visible between 4 and 5.4 km. Because the ice crystal content of the precipitation is large 
compared to the molecular backscatter, only small difference between rasv and inverted 
aerosol depolarizations is observed. The raw molecular depolarization shows small increase 
for the cloud, but after the inversion a constant 0.8% depolarization is observe^;. 


■t 
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The depolarization data for cirrus clouds obtained between August and November 
1993 have been analyzed. A summary of the observed depolarization ratios as a 
function of temperature are presented in Figure 30. In this study only clouds with 
scattering ratios greater than 0.5 were used. This was made in order to fully separate 
cloud depolarizations from noisy clear air depolarizations. The values of clear air 
depolarizations are low and affected by the photon counting statistics. Because the 
HSRL measurements have shown, that the scattering ratios of clear air aerosol and 
stratospheric aerosols can exceed values 1-3, the use of the scattering ratio of 0.5 does 
not guarantee a clear separation between clouds and clear air aerosols. On the other 
hand, the scattering ratios of the cirrus clouds can be below 1. Therefore, a visual 
separation between cirrus clouds and clear air aerosols is made and only the cloud 
altitudes are included to the study. In this study, 2 min averaging times for the data 
were used. By using a short averaging time, the errors due to temporal changes of 
the atmosphere were minimized. The atmospheric temperatures were obtained from 

: • n i i \ + 

radiosonde measurements and temperature intervals of 5 °C were used. Water clouds 
at cirrus cloud altitudes were separated from ice clouds based on the depolarization 
ratio values. Clouds with depolarization ratio values less than 15 % were classified as 
water clouds. 

The high depolarization values of the cirrus clouds are easily separable from the 
water cloud depolarizations. The average cirrus cloud depolarization varies from 
33% to 41% showing an increase towards the colder temperatures. Part of this can 
be expected to be from different shape, size, and orientation of the ice crystals at 
different temperatures 41 ’ 42 . The size and shape of the ice crystals haVe been found 
to be different at different temperatures and the crystals have been found to have 
a preferred orientation. A similar increase in cloud depolarization towards colder 
temperature was observed by Platt et a I. 43 . His study was made for midlatitude and 
tropical cirrus and the depolarization values for temperatures from-30°C'to -10 °C 
were lower, ranging from ~0.15 to ~ 0.25. The cirrus cloud depolarizations measured 
with the HSRL do not show any depolarization values below 20% for the teihperatures 
from -30 to -10 °C. The system used by Platt had a 2.5 mrad field of view, and 
the low depolarization values are possibly caused by multiple scattering from water 
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clouds. The HSRL measurements have shown that when the depolarization of a water 
cloud is measured with a system with ~ 1.0 mrad field of view, the effects of multiple 
scattering are large enough to cause depolarization values of 20%, gnd therefore the 
separation between water and ice becomes impossible. Also the HSRL. measurements 
show a substantial probability of observing water at temperatures from -30 to 0 °C. 

In 45% of the cases simultaneous observations of water cloud layers at the cirrus 
cloud altitudes were made. The Figure 30 shows, that supercooled water clouds 
at cirrus cloud altitudes are found at temperatures as low as -35 °C. Above 0°C 
temperature the low depolarization values indicate pure water, and the presence of 
cirrus disappears. The water cloud depolarizations are below 10 % and they contain 
the effects of multiple scattering. When the same clouds are studied with a NOAA- 
11 or NOAA-12 satellite, the simultaneous appearance of the water and ice cannot 
be noticed. On the satellite measurements, the cloud types are separated by using 
the information on the temperature. Therefore, if satellite data is used for cirrus 
optical depth studies, a supercooled water cloud layer mixed with ice cannot be easily 
separated and water will increase the optical depth value determined for the cirrus 
cloud. The depolarization ratio knowledge of the cloud measured with a lidar can be 
used to separate the water from ice and therefore separate optical depth measurements 
for both constituents can be performed. 

The initial measurements have shown that the HSRL is capable of measuring cirrus 
cloud particle sizes. These measurements together with phase function measurements 
can be used for further studies of the cloud particle size effects on the depolarization 
ratio 26 . > '■ ‘ 
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Figure 30, Cirrus cloud and supercooled water cloud depolarization as a function of temper- 
ature observed between August 2 and November 11, 1993. The depolarization values below 
10% are indication of water clouds and higher depolarization values indicate the presence 
of ice. 
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6.2 Measurements of scattering ratio, aerosol backscatter 
cross section, and optical depth 

Two examples from the measured scattering ratio, aerosol backscatter cross section 
and optical depth profiles for the November 11, 1993 dataset are presented in figures 
31 and 32. The cloud altitude and geometrical thickness of the clouds can be seen 
from the inverted aerosol profile. The inverted molecular signal shows the extinction 
in the signal through the atmosphere. This can be compared to the predicted Rayleigh 
return, calculated from the radiosonde data. The unsmoothed profiles are used. The 
Figure 31 shows a thin cirrus cloud case and the Figure 32 presents a water cloud 
together with a thin cirrus cloud. The water cloud is seen at 5.2 km altitude and the 
cirrus cloud can be found ranging from 7 to 9.2 km. The measured signals show the 
large dynamical range achieved in HSRL measurements: signal strength variations 
over four decades are seen. 

The scattering ratio is calculated from the ratio of inverted aerosol profile to the 
inverted molecular profile. The measured scattering ratios for the whole dataset vary 
from - 10~ 3 to -1000. The maximum scattering ratio is 10 for the case presented in 
Figure 31 and —300 for the water cloud case in Figure 32. 

The aerosol backscatter cross section is obtained from the scattering ratio and an 
estimate for the atmospheric density profile. The values of aerosol backscatter cross 
section change from — 10“ 6 m _1 sr“ 1 for cirrus and ice precipitation to — 10 _4 m _l sr“ 1 

for water cloud. The observed aerosol backscatter values of the clear air are between 
10~ 6 -and 10~ 7 m -1 sr -1 . 

The total optical depths are calculated from the ratio of the inverted molecular 
signal to the return predicted for the pure molecular scattering. They range from 0.5 
to 2.5 for this dataset. The Figure 31 shows a optical depth of 0.5 for the cirrus cloud 
and a optical depth 1.4 is observed for the case in Figure 32. The water cloud optical 
depth is 1 and the cirrus cloud gives an optical depth of 0.4. With the laser power 
used for this measurement, the optical depth of 2.5 is the upper limit of the current 
system. By increasing the laser output power larger optical depths can be measured. 
Because the cross talk between channels is small and it can be accurately corrected, 
and also because the 160 /irad field of view of the HSRL effectively suppresses multiple 

* l , , . r f ; « . 
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scattering, the optical depths inside a cloud can be measured. 

The accuracy of the optical depth measurements can be seen from the Figure 33. 
The leftmost graph shows the error in the optical depth measurement of a thin cirrus 
cloud and rightmost graph shows the error when a cirrus cloud of optical depth ~ 1 
is studied. The errors presented here are root mean square errors and they consist of 
photon counting, background correction, calibration, and wavelength tuning errors. 
The signals are spatially averaged with a 150 m filter and the filtering is taken into 
account in the error analysis. The accuracy of the optical depth measurements is 
limited by the photon counting statistics. This can be seen from the error bars, 
which are on the same order as the signal fluctuations. The accuracy of the HSRL 
measurements can be increased by increasing the signal strength and/or by using a 
constrained nonlinear fit to the inverted molecular profile. 
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Figure 31. The measured scattering ratio, aerosol backscatter cross section and optical 
depth profiles of a thin cirrus cloud between 6.5 and 10 km (November 117*1993 01:05-01:09 
UT). (a) The inverted aerosol and molecular profiles along with a calculated molecular 
profile. The difference between the calculated and measured profiles shows the extinction, 
(b) The scattering ratio values are ranging from ~0.1 (clear air between 3 and 6.5 km) to 
10 (cirrus cloud at 10 km), (c) The aerosol backscatter cross section of the cirrus cloud is 
~ 5-10 6 m *sr 1 and ~ for 10~ 7 m -1 sr -1 the clear air between altitudes 3 and 6.5 km. 
(d) The optical depth between 3 and 11 km is 0.55. The optical depth of the cirrus cloud 
is 0.4. 
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Figure 32. Scattering ratio, aerosol backscatter cross section and optical depth profiles of a 
thin cirrus cloud between 7 and 9.2 km and a water cloud with an ice crystal precipitation 
at 5.2 km (November 11, 1993 01:29-01:34 UT), (a) The inverted aerosol and molecular 
profiles along with the calculated molecular profile of the clear air. (b) Scattering ratio of 
the water cloud is ~ 300 and 10 for the cirrus cloud, (c) The observed aetosol backscatter 
cross sections are ranging from 10 _7 m“ 1 sr“ 1 for the clear air to 2-10 _4 m rl sr~ 1 for the 
water cloud, (d) The optical depth is 1 for the water cloud and 0.4 for the cirrus. 







74 




Figure 33. The root mean square errors in the optical depth measurement. A thin cirrus 
case (November 11, 1993 01:05-01:09 UT) (leftmost graph) and a cloud layer of optical depth 
of 1.4 (November 11, 1993 01:29-01:34 UT) (rightmost graph). The cases, are the same as 
in Figure 29. The error bars are of the same order as the signal fluctuations indicating that 
the optical depth measurements are mostly limited by the photon counting statistics. 
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6.3 Measurements of atmospheric temperature 

The temperature dependence of the Doppler-broadened molecular spectrum enables 
the measurements of the atmospheric temperature by a HSRL. The capability to 
measure the temperature would eliminate the possible errors due to the difference 
between the radiosonde reading and the current temperature over the lidar site. For 
this purpose, the temperature measurement capabilities of the University of Wisconsin 
HSRL were studied. The presence of clouds and strong layers of clear air aerosols will 
affect the shape of the molecular spectrum measured with the HSRL. The temperature 
measurement method presented in the following was used to test the capabilty of the 
HSRL to accurately measure the molecular spectrum width for layers with small 
aerosol content. 

Lidars have been used for the measurements of the atmospheric temperature pro- 
file by many lidar groups. The technique proposed by Strauch et aJ. 44 and Cooney 45 
allows calibrated temperature measurements by using the rotational Raman spectrum 
of nitrogen. With this technique about 1 °C temperature accuracy at low altitudes is 
achieved. Kalshoven et a/. 46 demonstrated a differential absorption lidar method for 
temperature measurements. They used 2 laser wavelengths and (^-absorption lines to 
measure atmospheric temperature up to 1 km altitude with 1 °C accuracy. Later En- 
demann and Byer 47 reported simultaneous measurements of atmospheric temperature 
and humidity with a continuously tunable IR-lidar. They used a three-wavelength 
differential absorption lidar technique and water vapor absorption lines. With this 
technique a 2.3 °C absolute accuracy was achieved. In addition to Raman and differ- 
ential absorption lidar techniques, Keckhut et aJ. 48 used Rayleigh scattering lidar to 
measure atmospheric temperature for altitudes 30-70 km. 

The temperature measurements made with a high resolution lidar have been re- 
ported by Alvarez et a l. 16 . Their temperature measurement is based on the two 
barium absorption filters with different bandpasses. Because the strength of the sig- 
nal received through an absorption cell is proportional to the width of the Doppler- 
broadened spectrum, the information of the signal strength together with a theoretical 
calculation for the Doppler- broadened Rayleigh-Brillouin spectrum can be used for 
the determination of the atmospheric temperature. Their latest measurements have 
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shown that only a 10 °C accuracy is achieved for profiles up to 5 km-. ■ 

The preliminary measurements of the atmospheric temperature made with the 
University of Wisconsin HSRL have been based on one iodine absorption filter. For 
temperature measurements the system transmission spectrum is measured by scan- 
ning the laser wavelength over the iodine absorption spectrum, similarly as in the 
system calibration scan. For one temperature profile, data from 5 calibration scans 
were averaged. This was done to increase the signal to noise ratio of the measurement. 
The measured profile was calculated by averaging the signal over a 300 m range with 
1 km steps. 

The signal from atmosphere and detected through the iodine absorption cell is a 
convolution of the Doppler- broadened molecular spectrum and the iodine absorption 
spectrum. The Brillouin modified approximation for the Doppler- broadened spec- 
trum was used to calculate the molecular line shapes at temperatures ranging from 
-70 to +30 °C with 1°C resolution. The calculated line shapes were convoluted with 
the measured iodine absorption spectrum. In order to define the atmospheric temper- 
ature at certain altitude, a least square fit was used to fit the measured profile to the 
calculated profile. The temperature that produced the best fit defined the tempera- 
ture of that altitude. Figure 34 shows an example of the received signal from 8 km 
altitude observed through the iodine absorption cell normalized by the signal observed 
with the channel without iodine absorption filter. The iodine absorption spectrum is 
shown as a reference. The modeled molecular profile is shown for the temperature 
that produced the best fit. The best fit was found at -47°C temperature. 
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Figure 34. The HSRL signal from 8 km altitude observed through the iodine absorption filter 
normalized by the signal measured simultaneously without the iodine absorption filter. The 
iodine absorption spectrum is plotted for reference. The modeled molecular transmission 
is shown for the temperature that produced the best fit between measured and calculated 
molecular transmissions. The best fit was obtained at -47°C temperature. 

The sensitivity of the molecular transmission of the iodine absorption filter to 
the width of the molecular spectrum is illustrated in Figure 35. The figure shows 
the effect of incorrect temperature to the fit. For this figure, the modeled molecular 
transmissions were subtracted from the modeled molecular transmission at -47 °C 
temperature. The temperature difference of 5 °C is displayed. The Figure 35 shows 
that a clear difference between temperatures is achieved, but because the differences 
are small, the accurate measurements of atmospheric temperature by the scanning 
technique are difficult to obtain. 
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Figure 35. The sensitivity of the molecular transmission to the temperature. The difference 
in modeled molecular transmission to the molecular transmission at -47°C is shown with 
5°C temperature steps. 

A temperature profile measured on February 27, 1994, in Madison between 00:00 
and 03:00 UT is presented in Figure 36. For the comparison the radiosonde tempera- 
ture profiles from nearest weather stations are presented. The stations at Green Bay 
(WI, 180 km northeast from Madison), Peoria (IL, 350 km south from Madison), and 
St. Cloud (MN, 450 northwest from Madison) provided a radiosonde profile at 00:00 
UT. The temperature values measured by the HSRL agree with the temperatures 
measured with the radiosondes. For the profile between 4 and 8 km the observed rms 
temperature differences are 2.97, 7.08, and 5.52 °C between HSRL and the weather 
stations. The rms difference between weather stations is 7.06 °C. For low altitudes, 
the largest difference between profiles is observed. This is expected because of the 
synoptic scale variations in weather conditions between different locations. 





79 


For altitudes between 6-8 km, a good agreement between HSRL measurement and 
the central Minnesota (St. Cloud) radiosonde profile is seen. These altitudes had a 
low aerosol content providing scattering ratio of ~ 0.02. The temperatures above 8 km 
show a big deviation from the radiosonde temparatures. This is due a strong aerosol 
layer, that disturbs the HSRL temperature measurement. The scattering ratio of the 
aerosol layer above 8 km was ~ 0.3. Also the measured temperature values for altitude 
between 2 and 3.5 km are colder than the radiosonde values. This is due to presence 
of a low level aerosol layer with scattering ratio of ~0.1 The presence of aerosols 
deepens the measured spectrum and therefore a fit into this spectrum underestimates 
the temperature. Therefore, if the temperatures are going to be measured in the 
presence of aerosols, the effect of the aerosol signal has to be separated from the 
molecular contribution. 
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Figure 36. A temperature profile obtained with the HSRL on February 27, 1994. The tem- 
perature profile obtained with radiosondes show the atmospheric temperature measured at 
the closest weather stations. A good agreement between HSRL and radiosonde observa- 
tions is observed between altitudes of 4 and 8 km. These are altitudes with a low aerosol 
backseat ter content. 
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7 Error analysis 

7.1 Error sources 

The signal detected with a photon counting system is a combination of backscatter 
signal, background light, and internal noise in the photodetector. The lidar equation 
(Eq. (1)) shows that the amount of received signal is determined by the number of 
transmitted photons, scattering cross section, backscatter phase function, area of the 
telescope, transmission of the receiver optics, and the sensitivity of the detectors. On 
the other hand, the HSRL measurements are dependent on the system calibration 
coefficients (Eqs. (17) and (18)). The errors in the HSRL measurements can be 
divided into 5 main categories: 

1. errors due to photon counting statistics 

2. errors due to changes in system calibration coefficients 

3. errors in the background subtraction 

4. errors due to a multiple scattering contribution 

5. errors in the atmospheric temperature profile used to compute the molecular 
density 


The error in the photon counting process is proportional to the square root of the 
measured signal. The errors due to photon counting statistics that are induced to the 
calibration coefficients were discussed in Chapter 5. 

The system detection efficiency is a combination of the system transmission and 
the photomultiplier quantum efficiency. The system detection efficiency affects the 
amount of detected signal and therefore, it is directly related to the photon counting 
statistics. A model to calculate the expected system transmission was developed. 
This model includes the spectral transmission of interference filters, etalons, and 
iodine absorption cell, transmissions and reflections of optical components (lenses, 
mirrors, cubes, etalon chamber windows) and photomultiplier quantum efficiencies. 
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The transmittance and reflectance values of different optical components were con- 
firmed with measurements. The calculated transmission efficiencies were compared 
with the measured values. For this study, a clear air case was chosen. Because the 
properties of molecular scattering are well known, the lidar equation for molecular 
scattering (Eq. 3) can be solved for the number of photons expected for the pure 
molecular return. 

In order to predict the expected molecular return, the calculated system trans- 
mission was convoluted with the calculated molecular spectrum at a 5 km reference 
altitude. The temperature for the Doppler- broadened molecular spectrum calculation 
was taken from a radiosonde profile. By calculating the ratio of measured to expected 
number of photons given by the lidar equation, the actual system transmission was 
estimated. The two way extinction between lidar and the backscatter volume was 
assumed to ~ 1. The comparison between expected and measured detection effi- 
ciency values is given in Table 2. This comparison was made without using the low 
resolution etalons. Each etalon decreases the tr ansmissi on by a factor of ~ 2. 


Table 2. h 

SRL detection efficiency 

Channel 

Calculated (%) 

Measured (%) 

Aerosol+Molecular 

0.78 

0.6 

Molecular 

0.44 

0.36 

WFOV 

4.4 

3.4 


The transmission efficiency values show ~ 25 % difference between calculated and 
measured values. A ~ 10% error can be made by assuming the extinction between 
lidar to 1. Also the uncertainty on the transmission and reflectance values of optical 
components affect the calculation. The largest uncertainty is the quantum efficiency 
of the photomultipliers. The value 0.12 is taken from the manufacturers specification, 
but only a 0.03 change in the quantum efficiency would correct the difference. The 
transmissions and reflectances of the receiver components that affect to the detection 
efficiency are listed in Table 3. The transmission and reflection numbers of some 
components are confirmed with a measurement and rest of the values are obtained 
from the manufacturers specifications. 
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Table 3. Transmission / Reflection losses in HSRL components. The source for the transmis- 
sion and reflection values is manufacturers specification (spec.) or a measurement (meas.). 


Channel 

Component 

R or T (%) 

Source 

Total (%) 

ALL 






Telescope mirrors (2) 

R = 90.0 

meas. 

R= 81.0 


Interference filter (1) 

T= 50.0 

meas. 

T=50.0 


Polarizing cubes (2/channel) 

T=99.0 

spec. 

T=99.0 



R=93.0 

spec. 

R=93.0 


PMT quantum efficiency 

QE= 12.0 

spec. 

QE= 12.0 

A+M, M 






Etalon chamber windows (3) 

T= 90.0 

meas. 

T= 73.0 


Mirrors (2 /channel) 

T= 90.0 

meas. 

T= 81.0 


Lenses (3 / channel) 

T=99.4 

spec. 

T=98.2 

A+M 






Beam splitter 

T=30.0 

meas. 

T=30.0 

M 






Beam splitter 

R=70.0 

meas. 

R=70.0 


Iodine absorption cell 

T (cell) =80.0 

meas. 

T=80.0 



T(mol.)~30.0 

calc. 

T=30.0 

WFOV 






Lens 

T=99.4 

spec. 

T = 99.4 

Total eff. 






Aerosol-1- Molecular Channel 



0.78 


Molecular Channel 



0.44 


WFOV Channel 



4.4 


The reflectivity of the telescope mirrors is limited by the reflectivity of fhe alu- 
minum coating. The poor transmission values of the etalon chamber windows and 
the mirrors have a big effect on the detection efficiency, because they lose ~ 41% of 
the received signal. These components are taken from the old system and. they are 
not optimized for the current operating wavelength. In order to increase the system 
performance, these components are going to be replaced in near future. 

The errors due to change in the system calibration coefficients are mainly caused 
by change in the transmission-receiver wavelength tuning and change in the system 
alignment. These errors can be calculated by a partial differentiation. 

A change in atmospheric temperature and pressure affects the amount Of detected 
molecular signal. Therefore, if the radiosonde profile is not valid for the measure- 
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ment conditions, an error is generated in the calculated molecular signal. The width 
of the Doppler-broadened spectrum directly affects the calibration coefficient C mm . 
This affects to the amount of molecular signal subtracted from the signal detected 
with the combined aerosol+molecular channel. In Brillouin spectrum and molecular 
backscatter cross section calculations, the effects of changing atmospheric .conditions 
are minimized by using the current radiosonde data from the nearest weather sta- 
tions. By comparing temperature and pressure readings from the closest radiosonde 
stations, an ~ 1-5 % difference between stations is observed. This due to the geo- 
graphic separation between weather stations. Because the lidar is located between 
these stations, the error is expected to be smaller. 

In addition to the atmospheric temperature and pressure changes, the effects of 
the wind on the measurement has to be considered because the high wind velocities 
can shift the Doppler- broadened spectrum with respect to the laser wavelength. A 
bulk shift in respect to the radiated spectrum is an indicator of wind, while the 
line broadening is associated with the random thermal motions of the molecules (and 
aerosols). If the line center of the Doppler-broadened spectrum is shifted in respect to 
the transmitted wavelength, the system calibration will be in error. The Doppler-shift 
for the backscatter signal can be calculated from 

A v v 

2 T = ~c' ( 24 ) 

where 

Ai/ = wavenumber shift 
2 v o — wavenumber of the laser line center 
v = molecular velocity 

c = velocity of light 


The wind velocities measured by a radiosonde can be used as an estimate for the 
horizontal wind. The maximum horizontal wind velocities are generally smaller than 
70 m/s at cirrus cloud altitudes, and the maximum vertical winds are normally below 
1 m/s in troposphere. Therefore, the maximum wind to the lidar pointing direction 
9 (9 — 4 degrees from the zenith) is vi = v x sin9 + v y cos9, which corresponds to a 
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velocity of 5.9 m/s. The shift in the spectrum is 19.8 fm, which is much smaller 
than the current 0.052 pm detection accuracy of the molecular spectrum. Also the 
horizontal wind speeds are generally much lower 49 , and therefore, the effects due to 
wind to the HSRL measurements are negligible, but caution is required if larger zenith 
angles of the pointing direction are considered. 

The uncertainty in transmitter-receiver peak transmission tuning can be seen as a 
shift in the transmitter wavelength in respect to the spectral transmission maximum 
of the receiver. This effects the amount of detected molecular signal, but if the drift is 
large, a change in the amount of detected aerosol is observed. The drift also affects the 
system transmission. In the current HSRL, the effects of a drift between transmitter 
and receiver wavelengths are minimized by locking the transmitter laser wavelength 
to the iodine absorption peak (see Chapter 5.2). The effects of the wavelength locking 
to the measured signal were also discussed in Chapter 5.2. 

In addition to the previous errors, the system calibration is affected by all mis- 
alignments of the system. The error due to a system misalignment can be divided 
into three different categories: 

1. alignment error in the receiver optics 

2. alignment error between transmitter and receiver pointing direction 

3. an error in the compensation for the beamsplitter polarization dependence 


A small alignment error (off focus, off axis or misalignment in system pointing 
direction) can lead situations, where part of the backscattered light is lost in the 
optics, detector, or atmosphere. This can easily lead into a range dependence of 
the received signal. This affects the inversion and furthermore, the measured optical 
depth. A misalignment also affects the total system transmittance and the signal to 
noise ratio of the measured signal. System pointing instabilities make the overlap 
correction of the received signal difficult. The overlap region is the close distance 
range where the overlap of the receiver field of view and the transmitted beam is 
incomplete. The effect of misalignment between tr ansmi tter end receiver pointing 
directions is largest for the overlap region. The effects vary daily depending from 
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the accuracy of the alignment between transmitter and receiver pointing direction. 
Large receiver field of views are less critical for the alignment errors. The signal 
measured with the large field of views of the WFOV channel can be used to align the 
small field of view channels. In the HSRL, the alignment between transmitter and 
receiver pointing directions is performed as follows. A corner cube is used to direct 
the outgoing laser beam to the receiver. The beam is aligned to the center of the 
receiver aperture by adjusting the turning mirror on the top of the telescope secondary 
mirror. Then the corner cube is removed and the signal from the atmosphere is used 
to verify the alignment. The signal of the large field of view of the WFOV-channel 
is compared to signal of the combined channel. If the channels are properly aligned 
and if the transmitter and receiver pointing directions are the same, the slopes of 
the signals are similar for the altitudes above the overlap region and therefore, the 
ratio between WFOV-channel signal to the spectrometer channel signal is constant. 
The alignment of the smallest WFOV channel aperture can be verified similarly by 
comparing the signal from the larger aperture to the signal from the smallest WFOV 
aperture. The effects of alignment errors on the system measurement accuracy are 
not included to the error analysis, because with a careful alignment the errors can be 
minimized. 

The accuracy of the background correction is mostly affected by the photon count- 
ing statistics. The effects of photon counting on the measured background are de- 
creased by summing at least 66 range bins together. The background correction 
method used in HSRL measurements was described earlier in Chapter 3.2.2. In ad- 
dition to the photon counting statistics, the accuracy of the background correction is 
affected by the slight tendency to overestimate the number of background counts for 
altitudes close to 33 km. Eventhough atmospheric density, and therefore the amount 
of Rayleigh scattering, decreases rapidly with the range, some of the signal that is 
measured into the upper range bins of the HSRL can still be due to Rayleigh scatter- 
ing. Therefore, when those range bins are used to calculate the number of background 
counts, the background is slightly overestimated due to the Rayleigh signal. 

The effects of the multiple scattering on the received lidar signal are usually ne- 
glected and the received signal is assumed to be single scattered. A lidar return from 
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water and thick cirrus clouds usually encounters large optical depths within a short 
distance of the cloud boundary and therefore, many of the received photons are likely 
to be from multiple scattering. The effect of multiple scattering to the lidar signal 
is dependent from the receiver field of view, cloud particle size, range from the lidar, 
and the optical depth of the cloud. The effect of multiple scattering can be seen in the 
lidar profile as a reduction in the rate of attenuation of the signal. In depolarization, 
the effect of multiple scattering will show up as an increase in the depolarization ratio 
towards the cloud top. The HSRL measurements have shown, that the effect of multi- 
ple scattering on the smallest field of view signals is small, because the depolarization 
ratio values observed for the water clouds are low. 

The current HSRL profiles are not corrected for multiple scattering contributions. 
An inversion program that accounts multiple scattering effects is currently in progress. 
The magnitude of the error caused by assuming the signal to be originated from single 
scattering depends from the field of view of the receiver, optical thickness of the cloud, 
size of the cloud particles, and penetration depth. The clear air aerosols and thin 
cirrus layers have a small or no multiple scattering effect. Because of the large particle 
size of ice crystals, the signal from thick cirrus clouds contains a multiple scattering 
contribution even for the small field of views. The multiple scattering from water 
cloud droplets is distributed to the larger field of views. The amount of multiple 
scattering on the measured signal can be estimated from the paper by Eloranta and 
Shipley 27 . 
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7.2 Effects of different errors on the measured optical pa- 
rameters 

The effects of different errors on the inverted aerosol and molecular return and to 
the measured optical parameters were calculated by partially differentiating Equa- 
tions (19), (20), and (4)-(13). The error analysis is made for the data obtained on 
November 11, 1993. A section of the data for a thin two-layer cirrus cloud structure 
is analyzed. The case is the same as shown in Figure 31. The thin cirrus cloud case 
is chosen for this study for two reasons. First, the backscatter signal from thin high 
altitude cirrus is small. Therefore, an error analysis made for the thin cirrus tests 
the measurement sensitivity of the system. Second, if accurate measurements of thin 
cirrus clouds can be made within short averaging time, then the HSRL can be con- 
sidered for studies of contrail formed cirrus. The short averaging time is important 
for the contrail studies, because they drift rapidly with the wind and only a limited 
number of samples from one contrail can be obtained. . < . v , . , . . . , . 

Information from cloud optical properties can be obtained by comparing the HSRL 
measurements with satellite observations. Thus, the averaging time of the HSRL 
data has to be close to the time resolution of a satellite image. The use of the 
short averaging times also prevents the smoothing of the the lidar signal so that the 
rapid changes in cloud structure in both time and altitude remain in the data. For 
this study, the 3 min averaging time was chosen to study the errors in the inverted 
aerosol and molecular returns, inverted aerosol and molecular depolarization ratios, 
backscatter ratio, optical depth, and aerosol backscatter cross section. Because the 
measurement of the backscatter phase function is difficult due to uncertainty in the 
extinction cross section determination, the 6 min averaging time was used for error 
analysis of the phase function measurement. 

The error analysis presented here shows the total errors together with the partial 
errors. All errors are calculated as mean square errors (see Equations 26-35) and 
presented as fractional root mean square errors (see Figures 37- 44). The partial 
errors in Figures 37- 44 present the effects of errors on the photon counting of the 
combined channel and the molecular channel, error in the molecular backscatter cross 
section per unit volume due to the uncertainty in temperature and pressure profiles, 
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and errors in the calibration coefficient determination. The error estimate calculations 
include the photon counting error, error in the background correction, error due to the 
uncertainty in the temperature and pressure profiles measured with a radiosonde, and 
errors due to the tuning of the laser wavelength. The errors due to system alignment 
and the error due to multiple scattering effects are not included to this error analysis. 
Also the error in the determination of the range R is negligible. 

The errors in background correction are assumed to be from photon counting 
statistics. The uncertainty of 2 % for the temperature profiles is used. The error 
in pressure profile is estimated to 1 mbar of the radiosonde pressure reading at each 
altitude. The error is due to the geographic separation between closest weather sta- 
tions. This is estimated from the radiosonde measurements from the closest weather 
stations. The uncertainty on the molecular spectrum calculation is estimated with 
the 2% uncertainty on the atmospheric temperature. The error on the calibration 
coefficient C mm determination is a combination of the photon counting error, un- 
certainty of the molecular spectrum, and the error caused by tuning of the laser 
wavelength. The accuracy of the calibration coefficient C am determination is limited 
by the photon counting statistics. The error in the molecular backscatter cross section 
per unit volume measurement depends on the errors in the atmospheric temperature 
and pressure. 

The error in the molecular profile can be given as 

(AN m (R))> = (A «.(*))» + + 


( 25 ) 
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which leads to equation 



Cam(fl) 2 AS a+m (fl) 2 

{3 m (R) — B m (R) — C arn (R) (S a+m (R) — 5 a+m (i?))) 2 

A S m (fl) 2 - 

(Sm(R) - Bm{R) - Cam{R) ( S a+m (R ) - B a+m (R))) 2 + 

A CarnW AC mm (fl) 2 

{CmmiR) “ Cgm(^)) (Cm m (J?) — C am (/?)) 2 

Cg m (i?) 2 Ag 0 + m (fl)2 

(^m (-R) B m (R ) — Cam(-R) (*5 a + m ( J R) — B a + m (R))) 2 

, 2fii 

(S m (fl) - - C^(fi) (S 0+ „,(fl) - B a+m (fi))) 2 W 


A5 m (i?) and AS’ a+m (i?) are the photon counting errors. A C am (R) and A C mm (R) 
the errors in the calibration coefficient determination. The background correction 
errors are given by A B m (R) and AB a+m (R). The error analysis made for the cali- 
bration coefficients show that C mm can be determined with better than 5 % accuracy 
and C ajn with a better than 2 % accuracy. 


The effects of different errors on the inverted molecular return are shown in Figure 
37. The errors in inverted molecular return are dominated by the photon counting 
statistics and the determination of the calibration coefficient C mm . The error due to 
the measurement accuracy of the aerosol signal is negligible. The errors in the C am 
determination have only a small effect on the inverted molecular return. 
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Figure 37. The effects of different errors on the inverted molecular profile (November 11, 
1993, 01:05-01:08 UT). Data from a thin cirrus cloud is analyzed by using a 3 min. averaging 
time. The measured molecular profile (rightmost graph) presents the signal variations as a 
function of altitude. The extinction due to the thin cirrus cloud is observed between 6.5 and 
10 km. The accuracy of the inverted molecular return determination is mostly limited by 
the photon counting statistics and the determination accuracy of the calibration coefficient 
C mm (leftmost graph). The large errors shown in the overlap region below 2 km are due 
to a focus error and incomplete overlap of the receiver field of view and transmitted laser 
beam. 
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The error in the measured aerosol profile can derived similarly to the molecular 
profile 


= ( l £§) 2 < AS ~ <*»* + ( S § j ) 2 + 

( l £ f § y ) 2 <**-<*»■ + (^ f 1 ) 2 < AC «"> 2 + 
(»)V^+(*)W 

This differential equation can be expanded to following form 


(27) 


(A NJJVV 
V N a(R) ) 


C mm {R) 2 AS a+tn (R) 2 

(S m (R) - B m (R) - C mm (R) ( S a+m (R ) - B a+m (R ))) 2 + 

ASJJRf 

(S m (R) - B m (R) - C mm (R ) ( S a+m (R ) - B a+m (R))) 2 + 

AC am (R) AC mm (R) 2 

(C am (/?) - C mm (R)) 2 (Cam(R) - C mm (R)) 2 

CmmjRfAB a+m (R ) 2 , 

(S m (R) - B m (R) - C mm (R) ( S a+m (R ) - B a+m {R ))) 2 

AB m (R ) 2 - 

( S m (R ) — B m (R) — C mm (R) ( S a +m{R ) — -B a+m (i ?))) 2 


(28) 


The accuracy of the inverted aerosol return is mostly affected by the determination 
accuracy of the inverted molecular return (photon counting statistics and C mm deter- 
mination accuracy) and the photon counting statistics of the combined aerosol and 
molecular channel (see Figure 38). For the cases of small aerosol backscatter content, 
large errors in the determination of the aerosol return are made when short averaging 
times are used. The errors are caused by subtracting a large amount of molecular 
signal from the combined channel signal that contains the strong molecular signal 
together with a small aerosol contribution. Therefore, the statistics of the molecular 
signal dominates the aerosol backscatter signal determination. The determination of 
clear air aerosols requires longer averaging times in order to achieve reliable results. 
On the other hand, the measurements of cloud aerosols can be done with ~ 1.0% 
accuracy (7-10 km). 



93 



Fractional Error Inverted Aerosol 


Figure 38. The effects of different errors on the inverted aerosol profile in case of a thin 
cirrus cloud (November 11, 1993, 01:05-01:08 UT). Averaging time of the data is 3 min. 
The measured aerosol profile (rightmost graph) shows the signal variation as function of 
altitude. The thin cirrus layer is observed between 6.5 and 10 km and a strong aerosol 
layer is seen between 0 and 3.5 km. The measurements of the aerosols are limi ted by the 
accuracy of the molecular profile measurements (leftmost graph). The 3 min averaging time 
provides 1-5% accuracy for thin cirrus cloud and strong aerosol layer measurements, but a 
longer averaging time is required for the measurements of the clear air aerosols. 
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Using the calculated errors for the measured aerosol and molecular profiles, the 
errors in the determination of the optical parameters (see Chapter 2) can be calculated 
as follows. 

The errors in the determination of inverted aerosol and molecular returns have a 
direct effect on the accuracy of the backscatter ratio (or scattering ratio) determina- 
tion (Figure 39). Therefore, the error in the backscatter ratio is 


(A SR{R)) 2 = 


( dSRW V 

\ dN a (R) ) 


(A N a (R)) 2 + 


f dSR(R) \ 2 

\dN m (R)) 


(AN m (R))\ 


(29) 


The effects of errors on photon counting, background correction, and calibration 
can be derived by combining the previous equation with the equations (26) and (28). 
Similar derivations of the differential errors can be made for the optical parameters 
given in following. 

For the cases of a low aerosol backscatter content, the errors in the backscat- 
ter ratio are dominated by the errors in the aerosol return determination: For a 
stronger aerosol backscatter return from a cloud, the errors due to molecular return 
determination are on the same order or higher than the errors due to aerosol return 
determination. The backscatter ratios of the thin cirrus and strong aerosol layers 
can be determined with better than 10% accuracy, but measurements of the clear air 
require longer averaging times. 
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Fractional Error Backscatter Ratio 


Figure 39. The effects of different errors on the backscatter ratio in case of a thin cirrus cloud 
(November 11, 1993, 01:05-01:08 UT). Averaging time of the data is 3 min. The backscatter 
ratio profile (rightmost graph) shows the ratio of the aerosol return to the molecular return 
as a function of altitude. Backscatter ratios ~0.1 to 20 are observed. , The errors in the 
backscatter ratio measurement (leftmost graph) are determined by the accuracy of the 
aerosol and molecular return measurements. For the altitudes with a low aerosol content, 
the error in backscatter ratio is limited by the accuracy of the aerosol backscatter return 
measurement. For the cirrus cloud the accuracy depends on the goodness of the molecular 
backscatter profile measurement. 
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The error in optical depth can be approximated as sum of error in the molecular 
scattering cross section per unit volume and error in the molecular return determina- 
tion. 

(Ar( *» 2 " 2 (S ) 2 + 2 GS ) 2 < A ^ m 

Error in the molecular scattering cross section per unit volume A is deter- 
mined by the uncertainties in the radiosonde based measurement of the atmospheric 
temperature and pressure. 

(A/W*)) 2 = (f^) W))’+ (“)V PWf (3D 

Error in optical depth measurement is dominated by the error in the determination 
of the molecular return and the uncertainty on the density profile measured by a 
radiosonde (see Figure 40). For this study, a 2 % error in the temperature profiles is 
assumed. The closest radiosonde measurement site is Green Bay (WI), which is ~ 
180 km northeast from the lidar. Because the weather conditions can vary between 
the lidar site and the closest weather station, larger errors in the temperature profile 
are possible. The effects of errors on the atmospheric density profile can be minimized 
by making radiosonde measurements on the lidar site. 

The figure shows, that with 3 min averaging time the cloud optical depths can be 
detected with ~ 10 % accuracy. This accuracy is sufficient when clouds with optical 
depths greater than 1 are measured. For situations where optical depth is less than 
1, a longer averaging time is required. 
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Figure 40. The effects of different errors to the optical depth in case of a thin cirrus cloud 
(November 11, 1993, 01:05-01:08 UT). Averaging time of the data is 3 min. The optical 
depth profile (rightmost graph) shows the variation of the optical depth as a function of 
altitude. Optical depth of 0.5 is measuered for the range from 3 to 12 km and optical depth 
of 0.4 is observed for the cirrus cloud between 6,5 and 10 km. The errors in the measurement 
of the optical depth below 6.5 km are dominated by the inaccuracy of the radiosonde profile 
(leftmost graph) and the photon counting statistics of the molecular channel. The error in 
the calibration coefficient C mm determination also has a significant effect on the total error. 
The optical depth of the thin cirrus cloud can be measured with ~ 10% accuracy. 
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The accuracy of the aerosol backscatter cross section measurement is limited by 
the accuracy of the molcular backscatter cross section per unit volume determination 
and the accuracy of the backscatter ratio. 



The effects of different errors on the aerosol backscatter cross section are shown 
in Figure 41. This figure shows, that the measurements are mostly limited by the 
photon counting statistics, but also the uncertainty on the fi m (R) determination has 
a significant effect. The aerosol backscatter cross sections of clouds and strong aerosol 
layers can be observed with better than 10 % accuracy, but the measurements of clear 
air aerosol backscatter cross sections require longer averaging time. 
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Fractional Error Backseat. Cross Sect. (1/(m*sr)) 


Figure 41. The effects of different errors to the aerosol backscatter cross section in case 
of a thin cirrus cloud (November 11, 1993, 01:05-01:08 UT). Averaging time of the data 
is 3 min. The aerosol backscatter profile is presented as a function of altitude (rightmost 
graph) and the backscatter cross section values range from 10 7 to 10 5 l/m. The aerosol 
backscatter cross section of the cirrus cloud (6-10 km) and the strong aerosol layer between 
1 and 3.7 km can be determined with 4-10 % accuracy, but the measurements of the aerosol 
backscatter cross section of the clear air require a longer averaging time (leftmost graph). 
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The error in the phase function is affected by the errors on determinations of the 
molecular scattering cross section per unit volume, the aerosol and molecular profile, 
and the extinction cross section. These errors can be further divided to the photon 
counting errors, errors in the calibration coefficient determination, and errors in the 
background subtraction. 


< 4tt ) 


t d EA 7r,-ft) \ 2 (d — t,R ) \ 

dfS^R) (&Pm{R)) 2 + qsrJr) 
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J \ 

f d P a (n,R)\ 2 


(A SR{R)) 2 + 


dPAR) 


mm* 


(33) 


Because the extinction section is a range derivative of the optical depth, the deter- 
mination accuracy of the molecular profile limits the phase function measurements. 
For this study, the accuracy of the phase function determination is estimated for a 
6 min section of the thin two-layer cirrus cloud. The statistics obtained within 3 min 
averaging time is not sufficient for the measurements of phase function profiles. 

The accuracy of the phase function value determination can be seen from Figure 
42. The cloud phase function can be observed with 10-20% accuracy when 6 min 
averaging time is used. By increasing the averaging time or the signal strength, 
accurate measurements of cloud phase function profiles can be made. 
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Fractional Error P( 180 )/ 4 *Pi 

Figure 42. The effects of different errors on the phase function in case of a thin cirrus cloud 
(November 11, 1993, 01:05-01:11 UT). The 6 min averaging time is used. The phase function 
profile is presented as a function of altitude and the average phase function of the cirrus cloud 
layer between 7,5 and 10 km is 0,02 (rightmost graph). The accuracy of the phase function 
measurements is determined by the photon counting statistics, determination accuracy of 
the calibration coefficients, and accuracy of the molecular scattering cross section per unit 
volume (leftmost graph). The accuracy achieved within 6 min averaging provides phase 
function measurements with ~ 20% accuracy for the cloud layer. 
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Error in the inverted depolarization ratio can be presented as a sum of errors in 
the parallel channel and the perpendicular channel signals. 


(A ^ W)2 - fe ) 2 + fe ) 2 <“> 


The Fig 43 shows that the accuracy of the depolarization measurements is mostly 
limited by the accuracy of the perpendicular channel signal. The errors in the perpen- 
dicular channel signal determination are dominated by the photon counting statistics. 
The Figure 43 shows that short averaging times provide accurate measurements of 
cloud depolarization, and therefore reliable separation between water and ice clouds 
can be based on the depolarization measurements of the HSRL. Also reliable depo- 
larization measurements of strong aerosol layers can be performed. 

Figure 44 presents the errors in the molecular depolarization ratio. The mea- 
surements of molecular depolarization ratio can be performed with better than 10 % 
accuracy for the altitudes between 0.8 and 4 km. Reliable measurements of molecular 
depolarization for higher altitudes require longer averaging times. By using long av- 
eraging times the effects of atmospheric temperature on the measured depolarization 
can be studied. 


i 1 *. . * 
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Fractional Error Aeros. Depolarization 


Figure 43, The effects of different errors to the inverted aerosol depolarization ratio in case 
of a thin cirrus cloud (November 11, 1993, 01:05-01:08 UT). Averaging time of the data is 3 
min. The depolarization profile shows the variations of the inverted aerosol depolarization 
as a function of altitude (rightmost graph). A ~ 40% cirrus cloud depolarization is observed 
(6.5 - 10 km) and the depolarization of the strong aerosol layer is ~ 5%. The measurements 
of the inverted aerosol depolarization ratio are limited by the accuracy of the perpendicular 
signal (leftmost graph). The depolarizations of clouds can be measured with ~1 % accuracy. 
The depolarizations of strong aerosol layers are obtained with better than 10% accuracy. 
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Figure 44. The effects of different errors on the inverted molecular depolarization ratio in 
case of a thin cirrus cloud (November 11, 1993, 01:05-01:08 UT). Averaging time of the data 
is 3 min. The depolarization profile shows the variations of the inverted molecular depolar- 
ization as a function of altitude (rightmost graph). A 0.8 % depolarization is observed from 
0.5 to 7 km. The measurements of the molecular depolarization are limited by the accuracy 
of the perpendicular signal measurement (leftmost graph). For altitude between 0.8 and 4 
km better than 10% accuracy is achieved. Measurements molecular depolarization for the 
higher altitudes require longer averaging times. 


- . • I 
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8 Conclusions 

The work presented in this thesis has shown that the installation of the iodine absorp- 
tion filter has substantially improved the performance of the University of Wisconsin- 
Madison High Spectral Resolution lidar (HSRL). The new system design that includes 
the iodine absorption filter, polarization, and multiple scattering measurement capa- 
bilities is shown. The high resolution etalon, that was used in the earlier system to 
separate between aerosol and molecular scattering, provided only 1:2 separation for 
the aerosol backscatter signal between channels. The iodine absorption filter has been 
shown to suppress the aerosol backscatter signal on molecular channel down to 0.08% 
and therefore the aerosol cross talk signal on the molecular signal can be easily re- 
moved even for optically thick clouds. In the etalon based system, the determination 
of the system calibration terms was sufficient to provide accurate inversions only for 
the clear air aerosols and thin cirrus clouds. The use of the iodine absorption filter has 
also improved the robustness of the HSRL, reduced the complexity of the system, and 
increased the optical transmission of the system, so that accurate measurements of 
the optical depth, backscatter cross section and phase function can be made. Also the 
simultaneous measurements of depolarization and multiple scattering are performed. 
The HSRL implementation shows a depolarization measurement technique that uses 
the one transmitter laser and one detector to measure both polarization components. 
Therefore, no calibration of the receiver is required. The multiple scattering mea- 
surements are realized with a separate channel that allows measurements of signal 
strength variations as function of field of view simultaneously with the measurements 
of the narrow field of view channels. 

The iodine absorption filter provides an absolute wavelength reference for the 
HSRL measurements. The iodine absorption line observed through a cell with 50% 
transmission on the line center is used for the wavelength locking of the HSRL trans- 
mitter laser. This provides stable operation over a long period of time without a 
need for frequent calibrations. Measurements have shown that the laser wavelength 
is maintained within ± 0.052 pm. 

The stability and reliability of the system have been tested by operating the io- 
dine absorption based HSRL at the University of Wisconsin-Madison campus. Start- 
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ing from July 1993, the HSRL has been routinely operated and data from different 
atmospheric conditions have been recorded. This dataset contains ~ 30 cirrus cloud 
cases. The measurements show that accurate measurements of optical properties of 
the atmosphere can be performed. The improved measurement accuracy has made 
possible to measure optical depth profiles inside the clouds. The current HSRL can 
be used to probe clouds that have optical depths up to ~ 3. This means that most 

of the cirrus cloud cases can be fully observed and the bases of thick water clouds 

1 ' » ■ ' • 

can be measured up to 300-500 m inside the cloud. This has been achieved by using 
the iodine absorption filter, high laser pulse repetition rate, small pulse energy per 
laser pulse, and very fast photon counting data system. The measurement accuracy 
of the HSRL is high enough to provide accurate measurements of optical parame- 
ters of clouds and strong clear air aerosol layers within 3 min averaging time, but 
the accurate measurements of clear air optical parameters require longer averaging 
times. The clear air optical parameters can be measured up to 35 km. The error 
analysis shows that accuracy of the HSRL measurements is mostly limited by the 
photon counting statistics. The system performance can be increased by increasing 
the system detection efficiency or/and increasing the transmitted laser power. The 
greatest improvement would be a photodetector with higher quantum efficiency and 
faster count rate capability. 

The depolarization data obtained by the HSRL shows the ability of thfe HSRL to 
distinguish between water and ice clouds. It is shown, that the 160/irad fi&d of view 
of the spectrometer channels effectively suppresses the multiple scattering effects on 
the measured depolarization ratio. Therefore, a reliable separation between water 
and ice clouds is possible. The HSRL measurements have shown, that traditional 
systems with 1-5 mrad field of views cannot reliably separate between water and ice, 
because the depolarization observed with a 1 mrad field of view is comparable to the 
ice depolarization within a small penetration depth from the water cloud base due 
to the multiple scattering. The error analysis shows, that the depolarization of the 
clear air aerosol layers can be observed with better than 10 % accuracy. The accuracy 
of the cloud depolarization measurements is better than 1 %. The measurements of 
depolarizations of weak aerosol layers and molecular backscatterer are limited by the 
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photon counting statistics. 

The study of the cloud depolarizations between August 2 and November 11, 1993 
shows that in 45 % of cirrus cloud cases simultaneous observations of supercooled 
water at cirrus cloud altitudes were made. The average cirrus cloud depolarization 
is shown to increase from 33% at -5°C temperature to 41% at -60°C. The observed 
behavior is different than observed by Platt et af. 43 . The depolarization observed by 
Platt et al. ranged from ~15% to ~40% between temperatures from -10 to -60 °C. 
The largest difference between the HSRL measurements and Platt’s measurements is 
observed for the temperature range from -30 to -10 °C, where Platt observed low cirrus 
cloud depolarizations. The low values of depolarization in Platt’s measurement are 
most probably caused by multiple scattering from supercooled water droplets, because 
the system used for this measurement had a 2.5 mrad field of view. The HSRL 
measurements show that supercooled water clouds have been found at temperatures 
as low as -40°C and pure water clouds have been found at temperatures above CPC. 
No water has been found at temperatures below -40°C and the presence of cirrus 
disappears at temperatures above 0°C. In the HSRL measurements, the cirrus cloud 
depolarization ratios for all temperatures are close to the values observed for the 
temperatures without supercooled water clouds. The small difference in the observed 
depolarization as a function of temperature may be a result of different shapes, sizes, 
and orientations of the ice crystals at different temperatures. 

The depolarization of the molecular backscatter is ~0.7-0.8%, when measured 
without the low resolution etalons. When the low resolution etalons are used, a 
~0.55-0.6% depolarization is measured. The depolarization measured for the signal 
from the Cabannes line and the rotational Raman lines is 1.5 % without, any spectral 
filters. The measured molecular depolarization value is larger than the expected 0.4 % 
depolarization of Cabannes line. The system filter bandpass admits a small fraction 
of the rotational Raman lines and blocks part of the Cabannes line and therefore 
an increase on the depolarization ratio is observed due to the presence of the highly 
depolarized rotational Raman lines. The model calculation for the depolarization 
transmission of the system show, that a 0.56-0.62 % depolarization is, expected for 
the case where no low resolution etalons were used. A 0.402 -0.425 % depolarization 
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is expected, when one or two etalon are used. The depolarization observed with 
the HSRL are larger than the expected values. The cause of the small difference 
between expected and measured depolarization values is currently unknown. The 
further analysis of the depolarization measurement accuracy of the HSRL requires an 
advanced study of effects of the iodine spectrum and rotational Raman lines to the 
depolarization. 

The HSRL measurements require a knowledge of the atmospheric temperature 
profile. The atmospheric temperature profile measured with the HSRL shows that 
the HSRL can be used to measure temperature with a high enough accuracy so that 
the measured temperature profiles can be used for the analysis of the HSRL data. 
Therefore, the requirement for radiosonde measurements of atmospheric temperature 
could be eliminated. Before the temperature measurements with the HSRL can be 
routinely performed, the effects of aerosol backscatter signal to the molecular trans- 
mission of the iodine absorption filter have to be removed. 

This study has provided an instrument basis for a design of a simple and robust 
lidar for the measurements of the optical properties of the atmosphere. The University 
of Wisconsin HSRL provides a unique instrument for the measurements of the cloud 
optical properties and the data measured with the HSRL provides useful information 
that can be used for the climate models that study the effects of clouds to the earths 
atmosphere. 
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